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Introduction:  Multiple  Myeloma  (MM)  is  an  aggressive  and  incurable  plasma  cell  malignancy 
often  characterized  by  IgH  EnhancerZ/WYC  (IgH/MYC)  translocations  that  drive  excess  levels  of 
the  c-Myc  oncoprotein.  Recently,  the  Bradner  laboratory  has  shown  that  inhibition  of  the  general 
transcriptional  co-activator  BRD4  with  a  selective  chemical  probe  (JQ1)  leads  to  dramatic  down 
regulation  of  c-Myc  expression  and  cell  death  in  MM  cell  lines  (Delmore  et  al.,  2011).  In  other 
tumors,  BRD4  inhibition  does  not  lead  to  down  regulation  of  the  lgH//WYC  translocation  gene, 
but  rather  causes  the  selective  down  regulation  of  other  key  cancer  genes  (Dawson  et  al.,  2011; 
Ott  et  al.,  2012;  Zuber  et  al.,  2011).  BRD4  is  a  BET  (bromodomain  and  extra-terminal)  family 
protein  (Filippakopoulos  et  al.,  2010)  that  binds  to  acetyl-lysine  residues  on  histones  and  other 
chromatin  associated  factors.  BRD4  is  a  key  co-activator  of  the  elongation  factor  P-TEFb  and 
has  been  shown  to  co-activate  transcription  through  co-operative  interactions  with  master 
regulator  transcription  factors  (Huang  et  al.,  2009).  P-TEFb  is  required  for  the  transcription 
elongation  of  essentially  all  active  genes  (Rahl  et  al.,  2010)  suggesting  a  general  role  for  BRD4 
in  broadly  co-activating  transcription.  Thus,  it  is  unexpected  based  on  current  paradigms  of 
mammalian  transcriptional  regulation  and  chromatin  structure,  how  inhibition  of  BRD4  can 
selectively  inhibit  the  transcriptional  activity  of  oncogenes  in  tumors.  Currently,  small  molecule 
inhibitors  of  BRD4  (including  derivatives  of  JQ1)  are  in  FDA  Phase  I  clinical  trials,  with  recent 
reports  of  positive  clinical  response  (AACR  2014).  As  such,  it  is  imperative  to  understand  how 
BRD4  regulates  gene  expression  and  the  mechanisms  by  which  inhibition  of  BRD4  leads  to 
dramatic  effects  at  a  small  subset  of  genes  in  a  cell  type  specific  manner.  This  study  is  focused 
on  understanding  mechanisms  of  JQ1  activity  in  MM  and  other  tumors  as  it  is  critical  for 
understanding  both  how  chromatin  regulators  function  in  cancer  and  how  general  transcriptional 
regulator  inhibitors  can  achieve  selective  effects  against  cancer. 

Keywords:  Chromatin,  Transcription,  Multiple  Myeloma,  MYC,  Therapeutics,  Gene  Regulation 
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Key  Research  Accomplishments  -  To  explore  the  mechanisms  by  which  inhibition  of  BRD4 
leads  to  selective  effects  on  oncogene  transcription,  we  have  undertaken  the  following  aims: 

Major  goals  and  accomplishments 

Aim  1.  To  map  Brd4  onto  the  transcriptional  and  epigenomic  landscape  of  MM 
Status:  Completed 

Reported  in  Loven  et  al.,  Cell  2013;  Chapuy  et  al.,  Cancer  Cell  2013,  Fulciniti  et  al., 
submitted 

Using  genome  wide  ChIP-Seq  (chromatin  immunoprecipitation  coupled  to  high  throughput 
sequencing)  approaches,  we  have  mapped  the  comprehensive  transcriptional  and  epigenomic 
landscape  of  MM  in  steady  state  and  in  response  to  treatment  with  increasing  doses  of  JQ1. 
Given  that  BRD4  was  present  at  all  active  promoters  and  enhancers,  it  was  unclear  how 
inhibition  of  such  a  broad  co-factor  could  yield  highly  specific  results.  In  several  JQ1  sensitive 
cell  lines,  we  integrated  global  BRD4  occupancy  data  with  comprehensive  chromatin 
landscapes  and  found  that  BRD4  was  disproportionately  found  at  only  a  handful  of  large 
enhancer  domains  that  spanned  tens  of  kilobases  and  were  also  highly  occupied  by  all  other 
chromatin  and  transcriptional  co-activators.  On  the  whole,  these  “super-enhancers”  made  up 
only  3%  of  the  cell’s  enhancers,  but  encompassed  40%  of  all  enhancer  bound  co-activators.  In 
each  tumor,  super  enhancers  associated  with  key  oncogenic  genes,  driving  the  not  only 
continued  cancer  cell  growth,  but  also  genes  specifying  tumor  cell  identity  itself. 

The  discovery  and  characterization  of  super-enhancers  first  published  in  this  work  (Loven  et  al., 
2013)  and  its  companion  paper  (Whyte  et  al.,  2013)  have  now  been  extended  and  further 
studied  in  more  than  40  publications.  Data  from  this  study  has  been  curated  and  made 
publically  available  on  the  Epigenome  gateway  browser  in  the  “Multiple  Myeloma  Epigenome 
Portal”  (http://epigenomegateway.wustl.edu/browser/).  Consequently,  we  have  developed  and 
made  openly  available  analysis  software  that  can  quantify  normalized  factor  occupancy  genome 
wide  to  identify  and  rank  super-enhancers,  genomic  regions  of  asymmetric  chromatin  co¬ 
activator  loading  (https://github.com/BradnerLab/pipeline).  These  methods  have  been  used  by 
more  than  100  laboratories  and  have  been  incorporated  into  online  databases  including 
dbSuper  (http://bioinfo.au.tsinghua.edu.cn/dbsuper/). 

In  subsequent  work,  we  have  employed  this  chromatin  mapping  framework  in  patient  MM 
samples  to  investigate  changes  to  the  enhancer  landscape  of  MM  induced  by  the  primary  tumor 
microenvironment  and  have  identified  candidate  primary  tumor  specific  super-enhancers.  We 
have  also  integrated  a  new  chromatin  mapping  technique  (ATAC-Seq)  that  identifies  putative 
cis-regulatory  elements  from  small  cell  samples  (Buenrostro  et  al.,  2013).  These  approaches 
have  been  extended  to  profile  transcription  factor  occupancy  in  the  MM. IS  cell  line  (Figure  1), 
and  also  to  profile  chromatin  landscapes  in  other  MM  cell  lines  (Fulciniti  et  al.,  in  preparation). 
Here  we  have  identified  distinct  regulatory  axes  governed  by  BRD4  and  the  cell  cycle  promoting 
transcription  factor  E2F  and  have  incorporated  this  theme  into  our  study  of  BRD4  cooperation 
with  other  transcriptional  regulators  (Aim  3).  Finally,  in  collaboration  with  others,  we  have 
extended  these  approaches  into  other  tumor  model  systems  including  Non  Hodgkin’s  lymphoma 
(Chapuy  et  al.,  2013),  triple  negative  breast  cancer  (Shu  et  al.,  submitted),  and 
medulloblastoma  (Lin  et  al.,  submitted). 

Aim  2.  To  examine  the  kinetic,  transcriptionai  response  to  BET  bromodomain  inhibition 
Status:  Compieted  and  extended  with  new  objectives 
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Reported  in  Loven  et  al.,  Cell  2013;  Chapuy  et  al.,  Cancer  Cell  2013;  Anand  et  al.,  Cell 
2013;  Brown  et  al.,  Molecular  Cell  2014;  Fulciniti  et  al.,  in  preparation 

From  data  generated  in  Aim  1,  we  hypothesized  that  asymmetric  occupancy  of  chromatin 
regulators  including  BRD4  at  oncogene  driving  super-enhancers  might  explain  their  sensitivity  to 
chromatin  co-activator  inhibition.  We  have  mapped  JQ1  induced  changes  in  gene  expression 
and  RNA  Pol  II  genomic  occupancy  in  both  a  time  and  concentration  dependent  manner  in  MM 
and  other  disease  model  systems  (Anand  et  al.,  2013;  Brown  et  al.,  2014;  Chapuy  et  al.,  2013; 
Loven  et  al.,  2013).  These  data  consistently  show  that  BET  bromodomain  inhibition  by  JQ1 
treatment  leads  to  a  global  decrease  in  transcriptional  activity,  specifically  a  decrease  in 
elongating  RNA  Pol  II.  This  inhibition  of  transcription  is  most  pronounced  at  super-enhancer 
proximal  target  genes  and  is  supported  at  the  chromatin  level  by  evidence  that  both  BRD4  and 
the  active  kinase  subunit  of  the  elongation  P-TEFb  are  preferentially  lost  at  super-enhancer  loci 
upon  JQ1  treatment.  In  MM,  JQ1  treatment  leads  to  rapid  downregulation  of  MYC  and  other 
super-enhancer  associated  oncogenes.  Interestingly  in  diffuse  large  B-cell  lymphoma  (DLBC), 
super-enhancers  associate  with  different  oncogenes  and  thus  JQ1  treatment  does  not 
preferentially  downregulate  MYC  transcription.  Instead,  in  DLBCL,  super-enhancers  associate 
with  other  prominent  B-cell  factors  including  BCL6  and  PAX5  that  are  strongly  downregulated 
upon  JQ1  treatment  (Chapuy  et  al.,  2013).  In  all  profiled  systems,  inhibition  of  chromatin  and 
transcriptional  regulators  including  BET  bromodomains  (BRD4),  the  BAF  complex  (BRG1)  and 
transcriptional  kinases  (CKD7)  caused  a  selective  inhibition  of  transcription  at  super-enhancer 
driven  oncogenes  as  compared  to  other  genes  (Chapuy  et  al.,  2013;  Chipumuro  et  al.,  2014; 
Kwiatkowski  et  al.,  2014;  Loven  et  al.,  2013;  Shi  et  al.,  2013). 

We  next  sought  to  explore  the  consequences  BRD4  inhibition  during  dynamic  cell  state 
transitions.  Previously,  we  had  quantified  BRD4  inhibition  response  in  static,  steadily  growing 
tumor  cell  models.  However,  MM  and  other  tumor  cell  states  are  often  maintained  through 
aberrant  oncogenic  signaling,  necessitating  an  understanding  of  the  role  of  BRD4  in  stimulus 
coupled  signal  transduction  response.  In  vivo,  MM  tumors  exists  in  a  pro-inflammatory 
environment  in  the  bone  marrow  that  is  supported  through  increased  NF-kB  signaling  and 
transcriptional  activity  (Demchenko  and  Kuehl,  2010).  NF-KB  transcription  factors  are  known  to 
directly  interact  with  BRD4  (Huang  et  al.,  2009)  and  we  hypothesized  that  activation  of  NF-kB 
transcriptional  response  subsequent  to  pro-inflammatory  stimulation  requires  BRD4  activity. 
Since  in  vitro  immortalized  MM  cell  line  models  fail  to  fully  articulate  in  vivo  inflammatory 
response,  we  utilized  the  well  characterized  NF-kB  responsive  primary  human  endothelial  cell 
system  to  investigate  basic  mechanisms  of  signal  dependent  BRD4  function. 

Treatment  of  endothelial  cells  with  pro-inflammatory  stimuli  caused  NF-KB  to  localize  to  the 
nucleus  leading  to  transcriptional  response  and  phenotypic  cell  state  transition  from  resting  to 
inflamed  and  activated  endothelium.  Surprisingly  we  found  that  NF-KB  activation  resulted  in  a 
rapid  and  comprehensive  remodeling  of  the  chromatin  co-activator  landscape,  with  more  than 
50%  of  all  BRD4  bound  super-enhancers  redistributing  to  sites  of  newly  acquired  NF-KB 
binding.  In  inflammatory  stimulated  endothelial  cells  co-treated  with  JQ1,  NF-kB  still  enters  the 
nucleus  and  binds  chromatin.  However,  BRD4  fails  to  redistribute  to  these  sites  of  NF-kB 
binding,  thus  blocking  super-enhancer  formation  and  consequently  abrogating  NF-kB  directed 
pro-inflammatory  transcriptional  response. 

These  data  suggest  a  crucial  role  for  BRD4  in  dynamic  cell  state  transitions  where  it 
communicates  signals  from  signaling  pathway  terminal  transcription  factors  to  RNA  Pol  II. 
Further,  JQ1  inhibited  transcription  of  NF-KB  directed  de  novo  super-enhancer  target  genes  an 
order  of  magnitude  more  than  genes  associated  with  resting  endothelial  super-enhancers  at 
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baseline.  This  observation  suggests  that  in  vivo  cell  populations  undergoing  dynamic  cell  state 
transitions  in  response  to  cell  signaling  are  most  likely  to  be  affected  by  JQ1.  Results  from  this 
study  of  BRD4  and  NF-kB  dynamics  were  published  in  Molecular  Cell  (Brown  et  al.,  2014)  and 
is  included  in  the  appendices  of  this  report. 

Aim  3.  To  explore  the  contribution  of  cooperative  binding  and  disproportionate  load  by 
BRD4  to  transcriptional  response 

Status:  50%  completed 

Reported  in  Fulciniti  et  al.,  in  preparation;  Lin  et  al.,  submitted 

The  initial  description  of  super-enhancers  in  MM  revealed  additional  features  of  super¬ 
enhancers  with  potential  utility  in  characterizing  tumor  epigenomes.  Super-enhancers  differ  from 
typical  enhancers  in  their  underlying  sequence  composition,  their  response  to  perturbation,  and 
their  ability  to  drive  high  levels  of  transcription  at  target  genes.  Towards  a  mechanistic 
understanding  of  these  features,  we  have  now  developed  experimental  and  computational 
approaches  to  map  transcription  factor  binding  sites  within  super-enhancer  loci  using  data 
derived  from  a  combination  of  H3K27ac  ChIP-Seq  to  map  super-enhancers  and  ATAC-Seq  to 
precisely  identify  high-resolution  transcription  factor  binding  sites.  From  this  data  we  are  able  to 
computationally  reconstruct  enhancer  gene  regulatory  networks  and  predict  transcription  factor 
regulated  gene  expression  programs.  Applying  these  methods  to  MM  cells,  we  observed  that 
many  cell  type  specific  transcription  factors  in  MM  were  regulated  by  a  super-enhancer  and  also 
bound  to  super-enhancers.  These  auto-regulatory  patterns  established  a  core  regulatory 
circuitry  in  which  a  small  number  of  transcription  factors  are  super-enhancer  associated,  and  in 
turn  bind  to  and  regulate  super-enhancers  (Figure  2).  Analysis  of  core  regulatory  circuitry  in  MM 
and  normal  cd19+  plasma  cells  revealed  many  shared  plasma  cell  transcription  factors  including 
XBP1  and  IRF  family  factors.  Differential  analysis  revealed  the  glucocorticoid  receptor 
transcription  factor  as  present  only  in  the  circuitry  of  MM  cells.  As  MM  are  primarily  treated  with 
dexamethasone  (a  glucocorticoid  receptor  agonist),  these  data  support  a  central  role  for 
glucocorticoid  receptor  signaling  in  driving  MM  oncogenesis  and  suggest  that  core  circuitry 
network  approaches  accurately  describe  tumor  cell  identity. 

We  have  also  used  methodology  developed  from  this  effort  to  elucidate  core  transcriptional 
regulatory  circuitry  in  other  poorly  classified  tumors.  In  collaboration  with  the  German  cancer 
research  institute  (DKFZ),  we  have  deployed  chromatin  enhancer  mapping  and  core  regulatory 
circuitry  inference  methods  to  elucidate  cellular  origins  of  different  medulloblastoma  subgroups. 
Medulloblastoma  is  clinically  classified  into  4  different  subgroups,  each  with  different  molecular 
pathologies.  However,  the  cellular  origins  of  individual  subgroups  are  poorly  understood.  Here 
the  use  of  core  regulatory  circuitry  inference  has  elucidated  the  developmental  and  cellular 
origins  of  the  previously  uncharacterized  Group  4  medulloblastoma.  A  manuscript  describing 
these  results  has  been  submitted  to  Nature  and  is  currently  in  revision  (Lin  et  al.,  2015 
submitted). 

Currently  we  are  extending  these  chromatin  profiling  and  computational  approaches  across  a 
large  number  of  hematopoietic  lineage  derived  normal  and  tumor  cell  populations.  To  date 
H3K27ac  ChIP-Seq  in  primary  patient  MM  samples  has  proven  technically  challenging. 
Although  ATAC-Seq  is  able  to  identify  transcription  factor  binding  sites  in  primary  samples,  it 
fails  to  adequately  map  super-enhancers,  necessitating  H3K27ac  mapping  approaches.  To 
overcome  technical  challenges  with  ChIP-Seq,  we  are  implementing  the  newly  developed  iChIP 
methodology  that  is  able  to  produce  ChIP-Seq  equivalent  data  from  small  (<  1  million  cell) 
patient  samples  (Lara-Astiaso  et  al.,  2014). 
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Finally,  analysis  of  the  transcription  factor  occupancy  in  relation  to  BRD4  in  MM  has  suggested 
new  therapeutic  strategies  to  target  the  disease.  Previously  we  reported  the  concomitant 
inhibition  of  MYC  and  E2F  activity  in  multiple  myeloma  (MM)  upon  treatment  with  the  BET 
bromodomain  inhibitor  JQ1  (Delmore  et  al.,  2011).  BET  bromodomains  (BETs)  are 
transcriptional  co-activators  that  occupy  active  promoters  and  enhancers,  but  are 
asymmetrically  localized  to  a  small  number  of  “super-enhancer”  domains.  JQ1  treatment  results 
in  disproportionate  displacement  of  BETs  from  super-enhancers  leading  to  potent  and  selective 
downregulation  of  super-enhancer  target  genes,  including  MYC  which  is  >  90%  depleted  after  6 
hours. 

In  contrast,  E2F  protein  levels  are  relatively  unperturbed  by  JQ1  treatment.  Instead,  JQ1 
treatment  inhibits  the  expression  of  super-enhancer  driven  upstream  regulators  of  E2F  including 
the  MYC  and  CCND2.  In  tumors  where  MYC  is  not  super-enhancer  associated,  JQ1  treatment 
fails  to  directly  downregulate  E2F  activity.  These  observations  suggest  that  the  MYC 
translocation  status  of  MM  plays  a  key  role  in  determining  the  sensitivity  of  E2F  activity  to  BET 
inhibition.  They  also  suggest  an  unexplored  collaboration  between  MYC,  E2F,  and  BETs  in 
maintenance  of  MM. 

In  both  MYC  translocated  and  non  MYC  translocated  MM,  we  mapped  the  global  occupancy 
patterns  of  E2F1  and  its  dimerization  partner  DPI.  We  next  utilized  chemical  and  genetic 
perturbations  to  investigate  the  functional  consequences  of  E2F  depletion  singly  or  in 
combination  with  JQ1  in  vitro  and  in  vivo.  Across  MM,  we  demonstrate  that  E2F  activity  is 
required  for  tumor  growth  both  in  vitro  and  in  vivo,  as  depletion  results  in  G1  arrest  (Figure 
3a, b).  We  find  surprisingly  that  DPI,  the  dimerization  partner  of  E2F1,  is  required  for  tumor 
growth,  and  that  DPI  expression  negatively  correlates  with  patient  outcome.  Global  chromatin 
analysis  reveals  distinct  regulatory  axes  for  E2F  and  BETs,  with  E2F  predominantly  localized  to 
active  gene  promoters  and  BETs  disproportionately  found  at  super-enhancers  (Figure  3  c-j).  As 
MYC  activates  E2F,  translocations  of  MYC  to  the  IgH  enhancer  place  both  super-enhancer  and 
E2F  driven  genes  under  BET  control.  Consequently  BET  inhibition  is  synergistic  with  E2F 
depletion  only  in  non-/g/-///WYC  translocated  MM  (Figure  3k). 

Our  results  implicate  E2F  as  a  dependency  in  MM  and  expose  a  vulnerability  to  BET  inhibition 
imparted  by  Igti/MYC  translocations.  In  non  IgH/MYC  translocated  MM,  E2F  inhibition  is 
synergistic  with  JQ1  treatment.  These  observations  suggest  targeting  of  E2F  as  promising 
therapeutic  strategy  in  MM.  A  manuscript  describing  these  results  is  currently  under  preparation 
for  2015  submission. 

Summary  of  reported  accomplishments:  During  this  reporting  period,  data  from  these  aims 
have  been  prepared  in  the  following  manuscripts. 

1)  Brown  et  al..  Molecular  Cell  2014 

2)  Lin  et  al.,  submitted  2015 

3)  Fulciniti  et  al.,  in  preparation 

4)  Wolf  et  al..  Trends  in  Cell  Biology  2014 

Opportunities  for  training  and  professional  development:  The  work  described  in  these  aims 
was  performed  in  collaboration  with  clinicians,  molecular  biologists,  and  computational 
biologists.  Through  collaborative  interactions,  I  have  gained  knowledge  and  experience  in  the 
clinical  management  of  cancer  and  translational  research  approaches.  Conversely,  I  have  been 
able  to  train  several  molecular  biologists  in  computational  bioinformatics. 
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Professional  development  during  this  reporting  period  include  the  presentation  of  invited 
seminars  at  the  Dana-Farber  Cancer  Institute,  Broad  Institute,  Memorial  Sloan  Kettering  Cancer 
Institute,  Baylor  College  of  Medicine,  and  the  National  Institutes  of  Health. 

Dissemination  of  resuits  to  communities  of  interest:  In  addition  to  the  aforementioned  public 
presentations,  work  from  these  aims  was  presented  to  disease  advocacy  groups  and  disease 
stakeholders  at  the  Dana-Farber  Cancer  Institute,  the  Leukemia  and  Lymphoma  Society,  and 
the  Multiple  Myeloma  Research  Foundation.  Emerging  themes  and  concepts  on  oncogenic 
transcription  regulation  by  the  transcription  factor  MYC  were  also  presented  in  a  peer  reviewed 
review  article  (Wolf  et  al.,  2014). 

Future  plans:  Research  priorities  for  the  next  reporting  period  include  submission  and  revision 
of  results  in  Fulciniti  et  al.,  and  the  preparation  of  a  manuscript  on  transcription  factor  networks 
in  MM  and  other  hematopoietic  malignancies. 
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Impact: 

The  discovery  and  characterization  of  super-enhancers  first  published  in  this  work  (Loven  et  al., 
2013)  and  its  companion  paper  (Whyte  et  al.,  2013)  have  now  been  extended  and  further 
studied  in  more  than  40  publications.  Data  from  this  study  has  been  curated  and  made 
publically  available  on  the  Epigenome  gateway  browser  in  the  “Multiple  Myeloma  Epigenome 
Portal”  (http://epigenomegatewav.wustl.edu/browser/).  Consequently,  we  have  developed  and 
made  openly  available  analysis  software  that  can  quantify  normalized  factor  occupancy  genome 
wide  to  identify  and  rank  super-enhancers,  genomic  regions  of  asymmetric  chromatin  co¬ 
activator  loading  (https://github.com/BradnerLab/pipeline).  These  methods  have  been  used  by 
more  than  100  laboratories  and  have  been  incorporated  into  online  databases  including 
dbSuper  (http://bioinfo.au.tsinghua.edu.cn/dbsuper/). 

The  discovery  of  BET  bromodomain  regulation  of  inflammatory  signaling  has  prompted 
additional  pre-clinical  research  on  therapeutic  strategies  to  target  chromatin  and  transcription 
regulation  in  vascular  disease  by  pharmaceutical  companies  developing  BET  bromodomain 
inhibitors  for  the  clinic  (Tensha  Therapeutics  &  Glaxosmithkiine). 

The  are  no  other  significant  impacts  to  report. 

Changes  and  Problems: 

We  have  experienced  technical  challenges  in  profiling  chromatin  in  primary  MM  patient 
samples.  This  is  largely  due  to  the  low  cellularity  of  MM  tumor  samples,  and  sample 
preservation  techniques  that  impede  chromatin  extraction.  To  overcome  these  problems,  we 
have  first  orthogonally  extended  our  approaches  to  other  tumor  systems  as  a  way  to  prototype 
new  analytical  methods.  We  have  also  extended  work  in  MM  cell  line  models  and  in  vivo 
xenograft  models  in  collaboration  with  other  MM  researchers.  Finally  we  are  in  the  process  of 
applying  new  low  cell  requirement  chromatin  profiling  techniques  including  the  newly  developed 
iChIP  methodology  that  is  able  to  produce  ChIP-Seq  equivalent  data  from  small  (<  1  million  cell) 
patient  samples  (Lara-Astiaso  et  al.,  2014). 

There  are  no  other  significant  changes  and  problems  to  report. 
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Products: 

Publications:  During  this  reporting  period,  data  from  these  aims  have  been  prepared  in  the 
following  manuscripts. 

1)  Brown  et  al.,  Molecular  Cell  2014 

2)  Lin  et  al.,  submitted  2015 

3)  Fulciniti  et  al.,  in  preparation 

4)  Wolf  et  al..  Trends  in  Cell  Biology  2014 

Presentations:  During  this  reporting  period,  this  work  was  presented  in  invited  seminars  at  the 
Dana-Farber  Cancer  Institute,  Broad  Institute,  Memorial  Sloan  Kettering  Cancer  Institute,  Baylor 
College  of  Medicine,  and  the  National  Institutes  of  Health. 

Conferences:  During  this  reporting  period,  this  work  was  presented  at  conference  poster 
proceedings  at  the  Keystone  Symposium  on  Epigenetics  and  Cancer  and  at  the  AACR  meeting 
“MYC:  From  biology  to  therapy”. 

There  are  no  other  significant  products  to  report. 
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Participants  and  Other  Coiiaborating  Organizations: 


Name: 

Jonathan  Brown,  M.D. 

Project  Role: 

Instructor 

Researcher  Identifier 

N/A 

Nearest  person  month 
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Contribution  to  Project: 

Dr.  Brown  contributed  experimental  analysis  and  direction  as  part  of  a 

collaboration  reported  in  Brown  et  al.,  2014. 

Funding  Support: 

Dr.  Brown  is  supported  by  the  NIH-K08  HL1 05678,  The  Watkins  Discovery 

Research  Award  and  The  Harris  Family  Award 

Name: 

Alexander  Federation,  Ph.D. 

Project  Role: 

Graduate  Student 

Researcher  Identifier 

N/A 

Nearest  person  month 

worked: 

4 

Contribution  to  Project: 

Dr.  Federation  provided  computational  analysis  as  part  of  a  collaboration 

reported  in  Aim  2  and  in  Lin  et  al..  Submitted  2015 

Funding  Support: 

Dr.  Federation  is  supported  by  a  Leukemia  and  Lymphoma  Society  SCOR, 

the  National  Science  Foundation,  and  NIH  grants  1R01  CA1 76745-01  and 

P01  CA1 09901 

Name: 

Mariateresa  Fulciniti  Ph.D. 

Project  Role: 

Instructor 

Researcher  Identifier 

N/A 

Nearest  person  month 

worked: 

12 

Contribution  to  Project: 

Dr.  Fulciniti  provided  experimental  analysis  as  part  of  a  collaboration 

reported  in  Aims  1  and  Aim  2. 

Funding  Support: 

Dr.  Fulciniti  is  supported  by  NIH  PO1-CA078378,  and  R01CA050947 

There  are  no  other  significant  participants  or  collaborating  organizations. 
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Appendices: 

1.  Figure  Legends 

Figure  1  Transcription  factor  and  chromatin  iandscape  of  muitipie  myeioma: 

a)  A  heatmap  showing  the  pairwise  similarities  between  spatial  occupancy  patterns  of  various 
transcription  factors,  chromatin  regulators,  and  chromatin  modifications  in  the  MM1.S  multiple 
myeloma  cell  line.  Individual  factors  are  hierarchically  clustered  to  group  spatially  similar  binding 
profiles. 

Figure  2  Transcription  reguiatory  circuitry  of  muitipie  myeioma: 

a)  Schematic  showing  methodology  of  constructing  transcription  factor  networks.  For  each 
super-enhancer  associated  transcription  factor  a  regulatory  IN  and  OUT  degree  are  calculated. 
For  any  transcription  factor  (TFg)  The  IN  degree  represents  the  number  of  other  super-enhancer 
associated  transcription  factors  that  binds  to  the  super-enhancer  of  TFg.  The  OUT  degree 
represents  the  number  of  super-enhancers  bound  by  the  transcription  factor  TFg. 

b)  Network  of  the  most  highly  interconnected  super-enhancer  associated  transcription  factors  in 
the  multiple  myeloma  cell  line  MM1.S.  Each  node  represents  a  super-enhancer  associated 
transcription  factor.  Each  edge  represents  binding  of  a  transcription  factor  to  the  super¬ 
enhancer  of  another  transcription  factor. 

c)  Comparative  network  analysis  of  various  hematopoietic  lineage  derived  cells.  Each 
comparison  shows  a  rank  ordered  bar  plot  of  the  largest  changes  in  total  regulatory  degree  (IN 
+  OUT)  for  each  cell  type  comparison.  At  each  transition,  transcription  factors  with  significant 
changes  in  total  regulatory  degree  are  highlighted. 

Figure  3  E2F  and  BET  bromodomains  establish  distinct  oncogenic  regulatory  axes  in 
multiple  myeloma: 

a)  A  panel  of  MM  cell  lines  were  infected  with  either  scrambled  (pLKO.1)  or  DP1 -targeted 
shRNA  and  selected  with  puromycin  for  72  hours.  DP1  mRNA  levels  and  cell  growth  were  then 
evaluated  by  qPCR  and  Thymidine  uptake  respectively.  The  results  are  presented  as  mRNA 
(gray  line)  or  cell  growth  (black  bars)  changes  from  cells  infected  with  pLKO.1.  Data  are  shown 
as  the  mean  values  ±  s.d.  of  triplicates. 

b)  In  vivo  evaluation  of  the  effects  of  DP1  knockdown  on  MM  cells.  Growth  curve  assess  tumor 
size  after  injection  of  MM.1S  cells  transduced  with  DP1-specific  shRNA  or  scrambled  control 
vectors  subcutaneously  into  the  right  posterior  flank  region  of  SCID  mice.  Data  are  shown  as 
the  mean  values  ±  s.d. 

c)  Gene  tracks  showing  RNA  Pol  II,  H3K4me3,  E2F1,  and  DP1  occupancy  at  the  E2F1  gene 
loci 

d)  Gene  tracks  showing  BRD4,  H3K27ac,  E2F1,  and  DP1  occupancy  at  the  BCL-xL  gene  loci. 
The  BCL-xL  intronic  super-enhancer  (SE)  is  annotated. 

e,f)  Pie  charts  showing  the  fraction  of  active  promoters  and  enhancers  bound  by  E2F1  and  DP1 

g)  Scatter  plot  showing  the  contribution  of  BRD4  enhancer  regulation  and  E2F  promoter 
regulation  for  all  active  genes  in  MM1.S.  The  y-axis  shows  enhancer  BRD4  signal  and  the  x- 
axis  shows  promoter  E2F  signal.  Units  are  in  reads  per  million. 

h)  Venn  diagram  showing  the  overlap  of  high  BRD4  SE  regulated  genes  and  high  E2F  promoter 
regulated  genes. 

i,j)  Functional  enrichment  categories  for  high  BRD4  SE  regulated  genes  (i)  and  high  E2F 
promoter  regulated  genes  (j) 

k)  Isobologram  of  anti-E2F  peptide  and  JQ1  treatment  effects  on  cell  proliferation. 
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Figure  1:  Transcription  factor  and  chromatin  landscape  of  multiple  myeloma 
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Figure  2:  Transcriptional  regulatory  circuitry  of  multiple  myeloma 
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Figure  3:  E2F  and  BET  bromodomains  establish  distinct 
oncogenic  regulatory  axes  in  multiple  myeloma 
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SUMMARY 

Proinflammatory  stimuli  elicit  rapid  transcriptional 
responses  via  transduced  signals  to  master  regula¬ 
tory  transcription  factors.  To  explore  the  role  of  chro¬ 
matin-dependent  signal  transduction  in  the  athero¬ 
genic  inflammatory  response,  we  characterized  the 
dynamics,  structure,  and  function  of  regulatory  ele¬ 
ments  in  the  activated  endothelial  cell  epigenome. 
Stimulation  with  tumor  necrosis  factor  alpha  pro¬ 
mpted  a  dramatic  and  rapid  global  redistribution  of 
chromatin  activators  to  massive  de  novo  clustered 
enhancer  domains.  Inflammatory  super  enhancers 
formed  by  nuclear  factor-kappa  B  accumulate  at 
the  expense  of  immediately  decommissioned,  basal 
endothelial  super  enhancers,  despite  persistent 
histone  hypieracetylation.  Mass  action  of  enhancer 
factor  redistribution  causes  momentous  swings  in 
transcriptional  initiation  and  elongation.  A  chemical 
genetic  approach  reveals  a  requirement  for  BET 
bromodomains  in  communicating  enhancer  remod¬ 
eling  to  RNA  Polymerase  II  and  orchestrating  the 
transition  to  the  inflammatory  cell  state,  demonstrated 
in  activated  endothelium  and  macrophages.  BET 
bronrKDdomain  inhibition  abrogates  super  enhancer- 
mediated  inflammatory  transcription,  atherogenic 
endothelial  responses,  and  atherosclerosis  in  vivo. 

INTRODUCTION 

Precise  control  of  inflammatton  is  essential  for  host  defense.  De¬ 
fense  against  pyogenic  infectbn  requires  rapid  activation  of  tis¬ 
sue  and  circulating  leukocytes  and  their  recruitment  by  activated 
endothelium.  However,  inflammation  is  also  integral  to  the  path¬ 
ophysiology  of  many  comrrxjn  and  life-threatening  illnesses. 
Acute,  high-grade  inflammation  accompanying  sepsis  features 


systemic  inflammatory  cell  activation  and  contributes  to  multi¬ 
system  organ  failure  and  death  (Medzhitov  et  al.,  2012).  Chronic, 
low-grade  inflammation  is  a  pathogenic  feature  of  autoimmune 
disorders  as  well  as  highly  prevalent  arxj  morbid  conditions 
such  as  diabetes  mellitus  arxj  atherosclerosis  (Libby  et  al., 
201 1).  As  such,  there  is  a  pressing  need  to  dissect  inflammatory 
signaling  for  the  elucidation  of  pathobgb  mechanisms  of  dis¬ 
ease  and  the  identification  of  targeted  therapeutb  interventions. 

In  inflammation,  a  primary  mode  of  bidirectional  cellular  com¬ 
munication  involves  one  set  of  cells  releasing  cytokines  to  acti¬ 
vate  surface  receptors  on  effector  cells.  Transduced  signals 
converge  on  activation  arb  translocation  of  inflammatory  tran¬ 
scription  factors  (Barnes  and  Karin,  1997).  A  central  pathway 
common  to  the  interaction  between  activated  teukocytes 
and  endothelium  is  tumor  necrosis  factor  alpha  (TNF-o^-medi- 
ated  signal  transductbn  to  nuclear  factor-kappa  B  (NF-kB)— a 
family  of  master  regulatory,  proinflammatory  transcription  fac¬ 
tors  canonically  defined  by  the  p50/p65  heterodimer  (Baltimore, 
2011).  Following  entry  into  the  nucteus,  NF-kB  birbs  to  DNA 
c/s-regulatory  elements  at  enhancers  and  prorrxiters,  prompting 
proinflammatory  transcription  pierce  et  al.,  1988). 

Gerxime-bound  nuclear  NF-kB  interacts  with  transcriptional 
coactivators  to  stimulate  transcription  at  multiple  steps  irxsluding 
the  remodeling  of  chromatin  as  well  as  the  initiation  and  elonga¬ 
tion  of  RNA  polymerase  II  (RNA  Pol  II;  Barboric  et  al.,  2001 ;  Kaik- 
konen  et  al.,  2013;  Natoli,  2(X)9).  NF-kB  recruits  and  interacts 
with  defined  chromatin  regulators  including  histone  acetyltrans- 
ferases  (P300),  histone  deacetylases,  and  epigenetic  reader  pro¬ 
teins,  such  as  BRD4  (Ashbumer  et  al.,  2001;  Huang  etal.,  2009; 
Zhong  et  al.,  2002).  Through  these  interactions,  NF-kB  engages 
in  crosstalk  with  chromatin  remodeling  machinery. 

BRD4  is  a  member  of  the  bromodomain  and  extraterminal 
domain  PET)  family  of  transcriptional  coactivators  and  elonga¬ 
tion  factors  (BRD2,  BRD3,  BRD4,  and  BRDT;  Dey  et  al.,  20(X); 
LeFtoy  et  al.,  2008).  At  active  genes,  BET  bronrxidomains  recruit 
the  positive  transcriptbn  ebngation  factor  complex  p-TEFb; 
Jang  et  al.,  2005;  Yang  et  al.,  2CX)5),  and  chromatin  remodeling 
factors  including  the  SWI/SNF  complex  (Shi  et  al.,  2013)  via  nx)- 
lecular  recognition  of  polyacetylated  histone  tails  (Mujtaba  et  al.. 
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2007).  Mechanistically,  TNF-a  or  lipopolysaccharide  (LPS)  stim¬ 
ulation  promotes  d  irect  acetylation  of  the  p65  subunit  (Lys31 0)  of 
NF-kB  by  P300  (Chen  et  al.,  2001),  promoting  adirect  interaction 
with  BRD4  through  twin  acetyl  lysine- recognizing  bromodo- 
mains  (Huang  et  al.,  2009).  This  interactbn  is  needed  for  produc¬ 
tive  NF-kB  transactivation  (Huang  et  al.,  2009),  suggesting  a 
central  role  for  BRD4  in  inflammatory  transcriptional  sigrtaling. 

Prior  research  from  our  group  and  others  has  identified  that 
BET  bromodomains  localize  genome-wide  to  promoter  arid  en¬ 
hancer  regions  (Anand  et  al.,  2013;  Chapuy  et  al.,  2013;  Loven 
et  al.,  2013;  Nicodeme  et  al.,  2010).  The  majority  of  enhancer- 
bound  BRD4  is  fourtd  within  a  small  number  of  massive  enhancer 
regions  termed  super  enhartcers  (SE).  Like  locus  control  regions 
or  stretch  enhancers,  SEs  concentrate  chromatin-bound  coacti¬ 
vators  to  gerres  essential  for  specialized  cellular  furtction  (i.e., 
immunoglobulin  production  in  plasma  cells),  and  lineage  specifi¬ 
cation  (i.e.,  germinal  center  differentiation;  Chapuy  et  al.,  2013; 
Loven  et  al.,2013;  Parker  etal.,  2013;  Whyteetal.,  2013).  Disrup¬ 
tion  of  SE  function  by  a  first  acetylated  lysine-competitive  small 
molecule  BET  bromodomain  inhibitor  from  our  group  known  as 
JQ1  (Filippakopoulos  et  al.,  2010)  suggests  the  mutability  of 
these  large  chromatin  structural  elements  (Chapuy  et  al., 
2013).  However,  the  role  of  SEs  in  the  control  of  dynamic  cell 
state  transitions  remains  unknown.  Recently,  BET  bromodomain 
inhibition  has  been  shown  to  abrogate  global,  maladaptive  tran¬ 
scriptional  programs  during  sepsis  and  heart  failure,  implicating 
BRD4  in  stress-induced  cell  state  transitions  (Anand  et  al.,  2013; 
Nicodeme  et  al.,  2010).  These  data  provided  a  rationale  for  inves¬ 
tigating  the  collaborative  roles  of  NF-kB  and  BRD4  in  regulating 
SEs  during  proinflammatory  activation. 

The  errdothelium  is  critical  to  the  initiation  and  propagation  of 
inflammatbn.  Endothelial  cells  (ECs)  prompt  leukocyte  recruit¬ 
ment,  adhesion,  and  trafficking  into  tissues,  thus  mediating  re¬ 
sponses  essential  for  many  inflammatory  disorders,  including 
atherogenesis,  in  which  activation  of  ECs  is  pathogenic  (Gim- 
broneetal.,  1990;  Ley  etal.,  2007).  Despite  these  important  roles 
in  disease,  global  studies  of  chromatin  structure  artd  function  in 
vascular  endothelium  have  to  date  not  been  undertaken.  In  this 
study,  we  investigate  the  role  of  BRD4  in  determining  the  inflam¬ 
matory  activation  of  ECs  through  NF-kB  and  SE  formation.  Here, 
we  provide  evidence  that  EC  activation  by  the  archetypal  proin¬ 
flammatory  stimulus  TNF-a  rapidly  depbys  NF-kB  to  enhancers 
and  promoters  genome-wide,  where  it  recruits  BRD4.  Through 
the  recruitment  of  BRD4,  NF-kB  establishes  new  SEs  coincWent 
with  the  surprising,  rapid  redistribution  of  BRD4  away  from  endo¬ 
thelial  resting  state  SEs.  Newly  established  NF-kB  SEs  are  prox¬ 
imal  to  arb  drive  cartonical  genes  of  the  inflammatory  response 
in  ECs,  including  key  effectors  of  chemotaxis,  adhesion,  migra¬ 
tion,  and  thrombosis.  BRD4  depletion  from  chromatin  through 
small  molecule  BET  bromodomain  inhibition  impedes  NF-kB- 
directed  SE  reorganization.  The  failure  to  form  proinflammatory 
SEs  preferentially  suppresses  SE-dependent  proinflammatory 
gene  transcription,  translating  into  functbnal  suppression  of 
key  TNF-a-induced  endothelial  responses  of  leukocyte  rolling, 
adhesion,  and  transmigratbn.  In  vivo,  we  find  that  BET  bromo¬ 
domain  inhibition  suppresses  atherogenesis— a  pathogenic 
process  predicated  on  inflammatory  endothelial  activatbn. 
Together,  these  data  establish  BET  bronrodomain-containing 


proteins  as  key  effectors  of  the  integrated  mammalian  inflamma¬ 
tory  response  through  their  rapid,  dynamic,  gbbal  reotganizatbn 
of  SEs  during  NF-kB  activation  and  suggest  SE  targeting  during 
inflammatory  cell  state  transitions  as  a  therapeutb  approach. 

RESULTS 

p65  and  BRD4  Establish  Super  Enhancers  during 
Proinflammatory  Stimulation 

To  explore  the  role  of  NF-kB,  BRD4,  and  SEs  in  the  acute  inflam¬ 
matory  activatbn  of  ECs,  we  activated  NF-kB  in  primary  human 
umbilical  vein  ECs  with  TNF-a,  a  canonical  prroinflammatory 
stimulus,  for  one  hour.  As  expected,  TNF-a  resulted  in  NF-kB 
nuclear  translocation  (Figures  1A  and  IB),  a  rapid  change  in 
EC  state  characterized  by  increased  monocyte  adhesion  (Fig¬ 
ure  1 C),  and  inductbn  of  adhesion  molecule  gene  expressbn 
including  E-selectin  {SELE)  and  vascular  cell  adhesion  molecule 
(l/CA/Wf ;  Figure  ID;  Ley  etal.,  2007).  The  recognized  capwcity  of 
BRD4  to  bind  acetylated  NF-kB  (Huang  et  al.,  2009),  suggests  a 
coactivator  role  for  the  BET  family  in  the  robust  p65- mediated  EC 
inflammatory  response  observed.  We  therefore  used  chromatin 
immunoprecipitatbn  coupled  with  high-thieughput  genome 
sequencing  (ChIP-seq)  to  define  p65  and  BRD4  genomic  occu- 
prarrcy  in  ECs  before  arb  following  acute  proinflammatory 
activatbn. 

InTNF-a-stimulated  ECs,  p65  enrichment  was  evident  at  pro¬ 
moters  (1 7.5%),  intragenic  (45.8%),  arb  intergenic  regulatory  se¬ 
quences  (36.7%;  Rgures  1 E  and  1 F).  Striking  colocalization  of 
BRD4  arb  p65  was  observed  by  gbbal  enrbhment  alignment 
arb  birbing  site  motif  analysis  (Figure  IF;  Figure  SI  A  available 
online).  TNF-a  treatment  prompts  dynamb  cobcalization  of 
p65  and  BRD4  to  enhancer  arb  promoter  regions  marked  by 
H3K27ac,  which  are  significantly  enriched  for  p65  consensus  se¬ 
quences  (Figure  S1A;  Matys  et  al.,  2006).  At  the  exemplary 
VCAMI  bcus,  TNF-a  stimulatbn  of  resting  ECs  for  one  hour 
increased  p>65  occupancy  at  both  prorrraters  and  uprstream 
enhancer  ebments  marked  by  acetylated  chromatin  (H3K27ac; 
Figure  1G).  Coirx:ident  with  these  events,  we  identified  recruit¬ 
ment  of  exceptbnally  high  levels  of  BRD4  at  discrete  hyperace- 
tylated  enhancer  elements  (Figures  IF  and  1G;  Figure  SIB), 
consistent  with  the  formatbn  of  de  novo  SEs  (SEs).  Focal 
BRD4  colocalizatbn  with  p65  was  observed  at  each  discrete 
peak,  with  complete  concordance.  Comparable  eviderrce  is  pro¬ 
vided  at  the  SELE  bcus,  where  TNF-a  stimulation  recruits  p65 
arb  high  levels  of  BRD4  to  a  gene  regulatory  region  completely 
devob  of  p65  and  BRD4,  augmenting  regional  hyperacetylatbn 
(Figure  SIC).  The  dramatic  remodeling  of  these  loci  in  one  hour  in 
TNF-a  stimulated  ECs  corroborates  the  robust  transcriptbnal 
activatbn  of  these  canonical  EC  inflammatory  gene  products 
(Rgure  1 D).  Notably,  typbal  enhancers  are  found  at  most  other 
EC  genes  as  exemplified  by  erbothelial  tyrosine  kinase  (7EK) 
and  serpin  peptbase  inhibitor  clade  HI  {SERPINH1),  where  the 
levels  of  p65  arb  BRD4  are  an  order  of  magnitude  lower 
compared  to  the  VCAM1  or  SELE  SE;  arb  TNF-a  does  not 
irbuce  mRNA  expressbn  (Rgure  1H;  Rgures  SID  and  SI  E). 

To  assess  the  genome-wide  distribution  of  SEs  during  the  EC 
inflammatory  cell  state  transition,  we  characterized  and  com¬ 
pared  the  enhancer  landscape  in  resting  arb  TNF-a-activated 
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Figure  1.  p65  and  BRD4  Genome  Binding 
during  Proinflammatory  Activation  in  ECs 

(A)  Images  of  ECs  ±  TNF  a  stained  for  p65  (red)  or 
DAPI  (blue)  (25  ng/ml,  1  hr). 

(B)  Western  blot  for  p65,  Ku  70,  and  tubulin  in 
cytosolic  (left)  and  nuclear  ^ight)  protein  fraction 
lysates  in  ECs  ±  TNF  a. 

(C)  Images  showing  adhesion  of  calcein  labeled 
THP  1  monocytes  to  ECs  ±  TNF  o  (25  ng/tnl,  3  hr). 

(D)  Bar  plots  showing  cell  count  normali2ed 
expression  levels  of  SELF  and  VCAM1  in  ECs  ± 
TNF  a  (25  ng/ml,  3  hr).  Enor  bars  represent  SEM. 
(^  Pie  chart  cf  p65  binding  site  dIstrBaution  in  EC 
genome  In  TNF  a(+). 

(F)  Heatmap  of  p65  (blue),  BRD4  ^ed),  and 
H3K27ac  (yellow)  levels  in  resting  ECs  and  after 
TNF  a  (25  ng/ml,  1  hr).  Each  row  shows  ±  5kb 
centered  on  p65  peak.  Rows  are  ordered  by  max 
p65  in  each  region.  ChIP  seq  signal  (rpm/bp)  is 
depicted  by  color  scaled  Intensities. 

(Q  and  H)  Gene  tracks  of  ChIP  seq  signal  for  p65, 
BRD4,  and  H3K27ac  at  the  VCAM1  arxf  TEK gene 
loci  In  untreated  (top)  or  TNF  a(+)  (bottom)  ECs. 
The  y  axis  shows  ChIP  seq  signal  (rpmArp).  The 
X  axis  depicts  genomic  position  with  TNF  a  gained 
typical  enhancers  (TE,  gray)  and  SEs  (SE,  red)  and 
promoter  regions  (whR^  marked. 

See  also  Figure  SI. 
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ECs  using  BRD4  ChIP-seq  data  sets.  When  ranked  by  increasing 
BRD4  enrichment,  347  and  271  SEs  were  identified  in  resting 
and  TNF-a-activated  ECs,  respectively.  These  SEs  comprised 
~7%  of  the  total  number  of  discrete  EC  enhancer  loci  (Figure  2A; 
Figure  S2A),  but  represented  more  than  one-fourth  of  the  total 
annount  of  enhancer  size  and  more  than  one-third  of  enhancer- 
bound  BRD4  (Figures  2A  and  2B).  Compared  to  typical  en¬ 
hancers,  SE  loci  are  significantly  larger  in  DNA  length,  total 
BRD4  signal,  and  signal  density  and  share  less  overlap  between 
resting  and  TNF-a-activated  ECs  (Figure  2C:  Figures  S2A  and 
S2B).  Following  TNF-a  stimulation,  the  absolute  change  in 
BRD4  total  signal  and  density  at  SEs  was  greater  compared  to 
typical  enhancers  (Figure  2D;  Figure  S2C).  We  observed  higher 


p65  total  binding  signal  and  density  at 
SE  loci  when  compared  to  either  typical 
enhancer  regions  or  active  gene  tran¬ 
scriptional  start  sites  (Figure  2E;  Fig¬ 
ure  S2D).  As  exemplified  by  the  VCAM1 
SE  locus  and  also  observed  globally, 
ECs  feature  dense  clustering  of  multiple 
regulatory  transcriptbn  factor  binding 
sites  known  to  be  involved  in  EC  proin¬ 
flammatory  responses  including  p65, 
p50,  ETS1  /2,  and  transcription  factor  3/4 
(TCF3/4:  Figure  2F:  Figure  S2E:  De  Val 
et  al.,  2008;  Masckauchan  et  al.,  2005). 
In  contrast,  typical  enhancer  sites  typified 
by  the  TEK  locus  possess  a  much 
lower  density  of  these  motifs  (Figure  2G: 
Figure  S2E). 

To  dissect  the  temporal  relattonship  between  p65  and  BRD4 
localization  to  enhancers,  we  next  performed  time-ranging 
chromatin  binding  studies.  ChIP  for  f>65  and  BRD4  foibwed 
by  real-time  PCR  centered  on  the  most  prominent  NF-kB  bind¬ 
ing  site  in  the  5'  VCAM1  and  CCL2  SEs  revealed  enrichment  of 
p65  five  minutes  after  TNF-a  stimulatbn,  with  peak  occupancy 
detected  by  30  min  (Figure  2H;  Figure  S2F).  BRD4  recruitment 
followed  the  identical  tem|X)ral  prattem  of  recmitment  at  these 
sites.  Inhibition  of  NF-kB  phosphorylatbn  arb  function  by  IkB 
kinase  inhibitbn  (BAY  11-7082,  “BAY”)  (Pierce  et  al.,  1997) 
completely  abrogated  TNFa-induced  p)65  and  BRD4  accumu¬ 
lation  at  both  NF-kB  sites  at  all  time  pxiints,  while  also  sup¬ 
pressing  VCAM1  gene  induction  (Figure  S2G).  As  expected. 
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Rgure  2.  p6S  and  BRD4  Establish  Super 
Enhancers  during  Proinflammatory 
Stimulation 

(A)  Ranked  plots  of  enhancers  defined  in  resting 
(top)  or  TNF «(+)  (bottom)  ECs  ranked  by 
increasing  BRD4  signai  (units:  ipm).  Enhancers  are 
defined  as  regions  of  BRD4  Chip  seq  binding  not 
contained  In  promotets.  The  cutoff  discriminating 
TEs  from  SEs  is  shown  as  a  dashed  line.  Genes 
associated  with  enhancers  that  are  considered 
typical  or  super  are  colored  gray  and  red, 
respectively. 

(B)  Pie  charts  displaying  characteristics  of  TE  and 
SE  regions  including  number  of  loci,  size,  and 
BRD4  signal. 

(C)  Boxplots  of  median  enhancer  length  (kb),  signal 
frpm),  and  density  (rpmAp)  in  TNF  a  gained  en 
hanceis.  Significanoe  of  the  difference  between 
distributions  determined  using  a  two  tailed  t  test. 
*”p<1  X  10 

(D)  Boxplot  of  absolute  change  In  BRD4  signal  in 
response  to  TNF  a  measured  at  all  enhancers  In 
TNF  a{  )  and  TNF  a(+).  Significance  of  the  differ 
ence  between  distributions  determined  using  a 
two  tailed  t  test.  *”p  <1x10  ’®. 

(E)  Boxplot  of  p65  binding  signal  (rpm)  at  all  active 
gene  promoters  (TSS),  TEs,  and  SEs  in  TNF 
a  treated  ECs.  Significance  of  the  difference  be 
tween  distributions  determined  using  a  two  tailed 
t  test.  **p  <  1  X  10  ®,  ***p  <1  X  10  ’®. 

(F  and  Q)  Schematic  of  transcription  factor  motif 
binding  sites  at  the  VCAM1  SE  ^ed  box;  F)  or  TEK 
TE  (gray  box;  G)  loci  In  ECs  treated  with  TNFo. 

(H)  Line  plots  of  kinetic  ChIP  PCR  showing 
enrichment  (percent  input  normalized  to  time  0)  of 
p65  and  BRD4  at  an  NF  xB  birxiing  site  In  the 
VCAM1  (lefQ  SE  and  TEK  TE  (right)  in  ECs  treated 
with  TNF  a  (25  ng/ml;  0,  5,  15,  30,  60  mkl).  The 
effect  of  cotreatment  with  vehicte  (top),  BAY  (NF 
kB  inhibitor,  mlddl^,  and  JQ1  (bottom)  Is  shown. 
See  also  Figure  S2. 
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cotreatment  of  ECs  with  JQ1  reduced  TNF-a-induced  enrich¬ 
ment  of  BRD4  at  these  NF-kB  binding  sites.  However,  in 
contrast  to  BAY,  JQ1  had  no  effect  on  immediate  early  p65 
recruitment  (Figure  2H;  Figure  S2F).  Compared  to  SEs,  typical 
enhancers  featured  less  TNF-a-induced  recruitment  of  both 
NF-kB  and  BRD4  (Figure  2H;  Figure  S2F).  Coilectiveiy,  these  re¬ 
suits  establish  that  p65  is  required  for  TNF-a-induced  recruit¬ 
ment  of  BRD4  to  SE  regtons  in  ECs. 


Tkiw (min)  NF-kB  Prowokes  Rapid  Global 

Redistribution  of  BRD4  at  Super 
Enhancers 

We  next  explored  the  dynamics  and 
funetbn  of  SEs  in  inflammation  in  both 
ECs  and  macrophages.  One  hour 
foliowing  TNF-a  stimulation,  we  ob- 
eo  served  the  formation  of  pronounced,  pro- 
totypicai  SEs  at  canonical  inflammatory 
genes,  such  cis  the  CCL2  chemokine 
figure  3A).  There,  a  density  of  upstream 
and  intragenic  enhancer  eiements  was  identified,  characterized 
by  regional  hyperacetylation  (H3K27ac)  and  peaks  of  BRD4 
enrichment  coinciding  with  focai  p65  binding  sites  (Rgure  3A). 
Functionaily,  rapid  SE  formation  was  associated  with  recruit¬ 
ment  of  RNA  Poi  II  (Figure  3A)  and  marked  transcriptionai  acti¬ 
vation  (Figure  S3A).  These  data  provide  a  demonstration  that 
the  canonicai  proinflammatory  stimuius  TNF-a  rapidly  induces 
de  novo  formation  of  SEs  estabiished  by  p65  at  proinflammatory 
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Figure  3.  NF  kB  Provokes  Rapid,  Global 
Redistribution  of  BRD4 
(A  and  B)  Gene  tracks  of  ChIP  seq  signal  ^pm/bp) 
for  p65,  BRD4,  H3K27ac,  H3K4me3,  and  RNA  Pol 
II  at  the  CCLZ  gene  (A)  or  SOX18  (B)  locus  in 
TNFa(  )  (top)  and  TNFo(+)  (bottom)  ECs. 

(C)  All  genomic  regions  containing  a  SE  in  TNFa(  ) 
and  TNF  a(+)  ECs  are  shown  ranked  by  logs 
change  In  BRDA  si^al  (treated  versus  untreated). 
The  X  axis  shows  the  log2  fold  change  In  BRD4 
signal.  Change  in  BRD4  levels  at  SEs  is  colored  by 
Intensity  of  change  fereen  to  red). 

(D)  Line  plot  showing  the  median  levels  of  p65 
binding  (rpm/bp)  at  SEs  In  either  TNF  a(  )  light 
blue  or  TNF  «(+)  dark  blue  conditions.  SEs  were 
ranked  by  change  in  BRD4  and  binned  ^0/bin). 
The  median  p65  level  was  calculated  in  each  bin. 
Error  bars  represent  95%  confidence  Intenrals  (Cl) 
of  the  median  determined  by  empirical  resampling. 
m  Horizontal  bar  plot  showing  the  ratio  of 
transcription  factor  motif  density  between  TNFa 
gained  and  TNF  a  lost  SEs.  Twenty  one  tran 
scriptlon  factors  are  displayed  whose  motifs  occur 
more  frequently  than  expected  based  on  dinucle 
otkfe  background  nrKxfel.  The  transcription  factor 
motifs  are  ranked  by  logs  fold  change  In  density 
between  TNF  a  gained  versus  TNF  a  lost  SEs. 
See  also  Figure  S3. 
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target  genes,  also  including  VCAMI  and  SELE  (Figures  1 G  and 
S1C). 

Unexpectedly,  the  enrichment  of  BRD4  at  inflammatory  EC 
SEs  was  associated  with  the  rapid,  reciprocal  depletion  of 
BRD4  at  resting  EC  SEs.  ECs  grown  in  culture  feature  347  ca¬ 
nonical  SEs,  including  SEs  associated  with  genes  critical  for 
noninflammatory  EC  function,  such  as  the  SOX18  transcription 
factor  implicated  in  vasculogenesis  (Matsui  et  al.,  2006).  In 
such  resting  ECs,  regtonal  enrichment  for  H3K27ac  is  observed 


at  the  SOX18  locus,  accompanied  by 
characteristic  BRD4  occupancy  and 
evident  enrichment  for  RNA  Pol  II 
throughout  the  gene  body.  In  the  inflam- 
matoiy  EC  state,  TNF-a  fails  to  drive 
NF-kB  to  the  30X18  locus,  RNA  Pol  II 
enrichment  is  markedly  diminished  and 
transcription  is  muted  (Figure  3B:  Fig¬ 
ure  S3A).  Only  1  hour  following  TNF-a 
stimulation,  BRD4  is  effectively  depleted, 
despite  persistent  hyperacetylation 
(H3K27ac)  of  the  30X18  SE  (Rgure  3B). 

To  explore  the  relevance  of  changes  in 
SEs  provoked  by  proinflammatory  activa¬ 
tion  across  the  geixrme,  a  systematic 
analysis  of  dynamic  alterations  in  SE 
formation  was  undertaken.  BRD4  en¬ 
richment  was  selected  here  as  a  pre¬ 
ferred  marker  for  SE  identity  owing  to  a 
concern  that  l-l3K27ac,  or  other  nrrare  bio¬ 
chemically  stable  enhancer  modifications 
(such  as  H3K4me1/2),  may  lag  behind 
BRD4  redistribution  globally  in  a  dynamic  cell  state  cf^nge,  as 
above  at  the  30X18  locus.  We  identified  a  dramatic  redistribu¬ 
tion  of  genomic  BRD4  occupancy  following  TNF-a  stimulation 
(Figure  3C).  Differential  enhancer  analysis  revealed  an  evident 
global  balance  in  chromatin-associated  BRD4,  but  importantly 
TNF-a  stimulation  resulted  in  the  reclassification  of  multiple 
enhancers  from  typical  enhancer  to  SE  (N  =  1 52)  and  from  SE 
to  typical  enhancer  (N  =  124;  Figure  3C).  Gains  in  BRD4  occu¬ 
pancy  were  strongly  and  directly  associated  with  site-specific 
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Figure  4.  NF  kB  Formed  Super  Enhancers 
Drive  Proinfiammatory  Transcription 
(A  and  B)  Bar  plot  of  average  change  In  elongating 
RNA  Pol  II  (A)  or  mRNA  expression  (B)  at  genes 
associated  with  TNF  a  gained,  TNF  a  lost,  or 
TNF  a  conserved  SEs  (ted),  TEs  (gray),  and  no 
enhancers  (black).  Error  bars  represent  95%  Cl  of 
the  mean  determined  by  empirical  resampling. 
Significance  of  the  differerxte  between  distrbu 
tions  determined  using  atwo  tailed  t  test,  "p  <  1  x 

10  ®,“’p<  1  X  10 

(C  and  D)  Change  in  elongating  RNA  Pol  II  in  the 
gene  body  region  of  genes  (C,  y  axis)  or  change  in 
mRNA  levels  (D,  y  axis)  are  plotted  ranked  by 
change  in  BRD4  at  proximal  SEs  (x  axis).  Dots 
represent  median  change  sarrpled  across  50 
evenly  distributed  bins  with  a  loess  fitted  line 
overlaid.  Change  In  BRD4  levels  at  proximal  SEs 
are  colored  by  intensity  of  change  (green  to  red). 
(E  and  F)  Metagene  representations  of  average 
RNA  Pol  II  Chip  seq  signal  (gray,  untreated:  black, 
TNF  a  treated)  in  units  of  rpm/bp  at  a  meta  com 
posite  of  target  genes  of  SEsthatare  gained  (E)  or 
lost  (F)  in  response  to  TNF  a  treatment.  Boxplots 
of  cell  count  normalized  expression  levels  are 
shown  to  the  right  of  each  metagene  in  arbitrary 
units  for  genes  with  associated  SEs  (gray,  un 
treated;  black,  TNF  a  treated)  that  ate  gained  (^or 
lost  (F)  In  response  to  TNF  a  treatment.  Signifi 
cance  of  the  difference  between  distributions 
determined  using  a  two  tailed  t  test.  **p  <  1  x 
10  ®.~*p<  1  X  10 

(Q  and  H)  Table  shewing  the  fumtlonal  categories 
of  selected  genes  that  are  targets  of  SEs  gained 
(Q)  or  lost  (H)  In  response  to  TNF  a  treatment. 

See  also  Figure  S4. 
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increases  in  p)65  binding  occupiancy  genome-wide  (Figure  3D). 
As  exfjected ,  TNF-a-gained  SEs  are  characterized  by  coordinate 
increases  in  BRD4  and  H3K27ac  enrichment  at  regions  of 
increased  p65  occupancy,  which  correiates  directiy  with  SE  for¬ 
mation  (Figures  3C  and  S3B).  Consistent  with  the  SOXfSreguia- 
tory  region,  lost  SEs  are  defined  by  a  significant  reduction  of 
BRD4  occupancy  foilowing  TNF-a  stimuiation  that  is  dispropor¬ 
tionate  to  minimai  changes  in  p65  and  modest  reductions  in 
H3K27ac  (Figures  3D  and  S3C).  Rnaliy,  the  p65  motif  was  found 
at  much  higher  density  in  TNF-a-gained  versus  TNF-a-lost  SEs 
regions  in  ECs,  further  suggesting  p65  direct  binding  to  clustered 


sites  in  TNF-a-gained  SEs  as  a  causai 
event  in  BRD4  SE  redistribution  during 
infiammatory  activation  (Figure  3E),  akin 
to  comparabie  “stretch  enhancers”  as 
described  by  Francis  Coliins  and  coi- 
ieagues  (Parker  et  ai.,  201 3). 

NF-KB-Formed  Super  Enhancers 
Drive  Proinfiammatory 
Transcription 

To  assess  the  consequences  of  SE  re¬ 
distribution  on  transcriptionai  output,  we 
integrated  genome-wide  ChIP-seq  data 
with  celi  count-normalized,  array-based  absoiute  gene  ex¬ 
pression  profiling  measurements  obtained  before  and  1  hour 
after  TNF-a  stimulation.  Compared  to  typicai  and  conserved 
enhancer  regions,  the  gain  or  ioss  of  SEs  provoked  by  TNF-a  re¬ 
suited  in  the  iargest  changes  to  RNA  Pol  II  occupancy  and 
expression  of  adjacent  genes  (Rgures  4A  and  4B).  This  globai 
reiationship  was  aiso  evident  when  specificaliy  comparing  the 
change  in  the  ievels  of  SE  constituent  BRD4  versus  change 
in  nearby  gene  eiongating  RNA  Poi  II  levels  and  expression  (Fig¬ 
ures  4C  and  4D).  Those  genes  positioned  near  TNF-a  gained 
SEs,  which  also  feature  BRD4  enrichment  at  their  promoters 
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suggestive  of  promoter-enhancer  communication,  demon¬ 
strated  marked  irtductionof  transcriptbnal  initiation,  ebngation, 
and  gene  expression  (Figure  4E;  Figure  S4A).  Functional  classF 
fication  of  TNF-a-gained  SE  marked  genes  reveals  known 
drivers  of  key  functional  facets  of  EC  inflammatory  responses: 
cytokine  signaling,  chemotaxis,  adhesion  and  migratbn,  arid 
thrombosis  (Figure  4G:  Figure  S4F).  In  contrast,  the  reciprocal 
loss  of  BRD4  at  resting  EC  SEs  resulted  in  a  proportbnate 
decrease  in  transcription  and  expression  of  nearby  genes  (Fig¬ 
ures  4A,  4B,  and  4F),  such  as  THBD  and  genes  involved  in  angio¬ 
genesis  and  endothelial  barrier  function  (Figure  4H;  Figure  S4G). 
Notably,  the  cohort  of  lost  SE  genes  featured  marked  decreases 
in  transcriptional  initiation  and  elongation  by  RNA  Pol  II  enrich¬ 
ment  metagene  analysis  in  advance  of  changes  in  prorrxjter 
modification  (H3K4me3  enrichment:  Rgure  S4B).  Genes  with 
conserved  enhancers  show  minimal  change  in  expressbn  and 
serve  prathways  that  govern  homeostatb  furxitbn  in  ECs  (Figures 
4A  and  4B;  Figure  S4H).  Globally,  this  chromatin  restructuring 
results  in  a  strong  inductbn  of  proinflammatory  SE  driven  tran¬ 
scription  compared  to  typical  enhancer  associated  genes.  The 
62  TNF-a-gained,  SE-associated  genes  comprise  only  ~8%  of 
all  genes  with  >2-fold  increase  in  mRNA  expressbn,  but  account 
for  ~60%  of  the  total  increase  in  upregulated  gene  expression 
and  ~20%  of  the  increase  in  cellular  mRNA  within  3  hr  of  TNF- 
a  stimulation  (Figures  S4C-S4E).  These  data  provide  discrete 
examples  and  global  evbence  of  dynamb,  functbnal  remodeF 
ing  of  enhancer  factors  even  preceding  the  structural  decommis¬ 
sioning  of  abandoned  enhancers. 

NF-xB-Formed  Super  Enhancers  Drive  Proinflammatory 
Gene  Expression  in  a  BET  Bromodomain-Dependent 
Manner 

Disruption  of  SEs  by  disrupting  enhancer  factors,  such  as  BRD4, 
has  been  observed  to  selectively  influence  the  expression  of 
genes  associated  with  SEs  (Loven  et  al.,  2013).  In  cancer,  we 
have  observed  that  competitive  displacement  of  BET  bromodo- 
mains  from  nuclear  chromatin  provokes  coordinated  inhibition  of 
the  MYC  transcriptional  program,  often  associated  with  downre- 
gulation  of  MYC  itself  (Chapuy  et  al.,  2013;  Delmone  et  al.,  201 1 ; 
Zuber  et  al.,  201 1).  Using  a  chemical  genetb  approach,  we  as¬ 
sessed  the  role  of  BET  bromodomains  in  the  rapb  transcriptional 
response  of  SE-associated,  proinflammatory  genes  in  TNF- 
a-stimulated  ECs.  Small  molecub  probes,  such  as  the  BET  bro- 
modomain  inhibitor  JQ1 ,  are  partbulariy  appealing  in  the  study 
of  dynamb  processes,  because  they  offer  precise  temporal 
perturbation  of  the  biological  system  (Frye,  2010;  Strausberg 
and  Schreiber,  2003).  We  therefore  performed  a  dynamic, 
genome-wbe  analysis  of  BET  bromodomain  inhibitbn  on  inflam¬ 
matory  SE  integrity,  global  chromatin  structure,  gene  expres¬ 
sion,  and  EC  postinflammatoiy  functbn. 

Rrst,  we  characterized  the  effect  of  JQ1  on  EC  chromatin 
structure  and  function  immediately  foibwing  TNF-a  exposure. 
ECs  were  treated  with  JQ1  (500  nM)  to  dispiace  BET  bromodo¬ 
mains,  then  stimulated  with  TNFa  for  1  hour.  Chromatin  from 
treated  and  untreated  ECs  was  subjected  to  ChIP-seq  for  pro¬ 
moters  (H3K4me3),  enhancers  (H3K27ac),  RNA  Pol  II,  the  p65 
transcriptbn  factor,  and  the  BRD4  coactivator.  At  the  SELE  lo¬ 
cus,  BET  inhibitbn  had  no  effect  on  TNF-a-mediated  recruit¬ 


ment  of  p65  (Figures  5A  and  5B).  Flowever,  vX31  depteted 
BRD4  resulting  in  abrogated  SELE  expression  assessed  by 
decreased  RNA  Pol  II  enrichment  (Figure  5A  and  cell  surface 
protein  levels  assessed  by  flow  cytometry  (Figure  S5E).  Com- 
p>arable  observations  are  evident  at  the  CCL2  and  IRAK2  loci, 
where  BET  inhibition  selectively  displaces  BRD4  leading  to 
impaired  transcription  induction  and  eiongatbn  by  RNA  Pol  II 
(Figures  S5A  and  S5B). 

Gbbal  analysis  of  TNF-a-stimulated  ECs  revealed  preferential 
loss  of  BRD4  at  TNF-a-gained  SEs  compared  to  typical  en¬ 
hancers,  with  minimal  effect  on  TNFa-induced  f>65  birbing  or 
l-l3K27ac  levels  (Figure  5B).  TNF-a-gained  SEs  possessed  a 
diminished  capacity  to  drive  proinflammatory  transcription  initi¬ 
ation  and  elongation  in  the  presence  of  the  BRD4  inhibitor, 
JQ1  (Figure  5C),  and  these  dynamb  effects  on  RNA  Pol  II  enrich¬ 
ment  were  independent  of  alterations  in  promoter  modification 
by  differential  metagene  analysis  of  FI3K4me3  enrichment  at 
the  transcriptional  start  sites  (Figure  S5C).  Functionaliy,  prefer¬ 
ential  JQ1  -induced  bss  of  transcription  at  TNF-a-gained  SE- 
associated  genes  culminated  in  more  potent  suppression  of  their 
proinflammatory  gene  expression  program  relative  to  genes 
driven  by  typical  enhancers  (Figure  5C:  Figure  S5D).  Treatment 
of  ECs  with  TNF-a  and  JQ1  suppressed  the  maximal  mRNA  in¬ 
ductbn  of  SE-associated  genes  (FS3,  CCL2,  VCAM1)  at  a  lower 
concentration  of  JQ1  and  to  a  greater  degree  compared  to 
TE-associated  genes  (LOX,  TEK,  NURP1\  Figure  5D).  Together, 
these  data  support  a  nxidel  in  which  BET  bromodomains 
mediate  dynamic  and  immediate  inflammatory  EC  response 
gene  transcriptbn,  by  facilitating  chromatin-dependent  signal 
transductbn  from  NF-kB  to  RNA  Pol  II. 

To  determine  whether  JQ1  transcriptionai  effects  resulted 
from  specific  engagement  of  BETs  at  nuciear  chromatin,  we 
spatially  localized  the  JQ1  molecule  genome-wide  using  a  new 
biotechnology  called  Chem-seq  (Anders  etal.,  2014).  This  tech¬ 
nique  maprs  the  interactbns  of  small  molecules  with  chronnatin  in 
the  human  gerK>me,  in  this  instance  by  using  a  retrbvable  syn¬ 
thetic  derivative  of  JQ1  (biotin-JQ1).  Notabiy,  despite  high  levels 
of  acetylated  HK27ac  at  enhancers,  we  detected  no  bbtin-JQI 
occupancy  \«ith  Chem-seq  at  SELE  and  other  SEs  (Figure  SE; 
Figures  S5F  and  S5G)  in  resting  ECs.  Rather,  at  these  loci, 
biotin^QI  birbing  to  chromatin  perfectly  cobcalized  with 
BRD4  spatially  arb  temporally  following  TNF-a-stimulation. 
Genome-wide  analysis  demonstrated  a  dose-dependent  reia- 
tionship  between  biotin-JQ1  localization  arb  the  decrease  in 
genomic  BRD4  occupancy,  not  H3K27ac,  in  ECs  treated  with 
TNF-a  +  JQ1  compared  to  TNF-a  alone  (Rgures  5F  and  5G). 
Taken  together,  these  data  reveal  that  JQ1  directly  targets 
BRD4  during  proinflammatory  activation  in  ECs.  Preferential 
loss  of  BRD4  at  inflammatory  SEs  by  BET  bromodomain  inhibF 
tion  serves  mechanistically  to  underscore  the  observed  selective 
effects  on  transcription. 

To  assess  the  generalizabiiity  of  these  observations,  we  tested 
whether  the  BRD4-deperbent  formation  of  proinfiammatory  SEs 
is  retevant  for  other  immune  effector  cells.  We  integrated  robust, 
publically  available  acetyFhistone  data  (H4K12ac)  with  RNA  Pol 
II  ChIP-seq  data  in  LPS-stimulated  macrophages  (Nbodeme 
et  al.,  2010),  to  map  SEs  and  explore  the  effect  of  BET  inhibition 
on  enhancer-mediated  transcriptional  signaling.  FI4K12ac  is  a 
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Rgure  5.  NF  kB  Formed  Super  Enhancers 
Drive  Proinflammatory  Gene  Expression  in 
a  BET  Bromodomain  Dependent  Manner 
(A)  Gene  tracks  of  ChIP  seq  signal  (rpmA>p)  for 
p65,  BRD4,  H3K27ac,  H3K4me3,  and  RNA  Pol  II  at 
the  SELB  locus  in  TNF  a  treated  ECs  ootreated 
with  vehicle  (top)  or  JQ1  (bottom). 

P)  The  mean  log:  fold  change  in  H3K27ac  (yellow), 
BRD4  (ted),  and  p65  (blue)  Chip  seq  signal  in  TNF 
a  treated  cells  ±JQ1  at  either  TEs  or  SEs  gained  In 
response  to  TNF  a  treatment.  Error  bars  represent 
95%  Cl  of  the  mean  determined  by  empirical  re 
sampling. 

(C)  Metagene  representations  of  average  RNA  Pol 
II  Chip  seq  signal  (black,  TNF  a  treated;  red,  JQ1 
treated)  in  units  of  rpmAip  at  a  meta  composite  of 
target  genes  of  SEs  gained  In  response  to  TNFa 
treatment.  Boxplots  (right)  show  cell  count 
normalized  expression  levels  In  TNF  a  (25  ng/ml, 
3  hr)  treated  ECs  ±  JQ1.  SIgnificanoe  of  the  dif 
ference  between  distributions  determined  using  a 
two  tailed  t  test.  ***p  <1x10 

(D)  Line  plots  of  mRNA  levels  (qRT  PCR)  of  three 
representative  genes  associated  with  TEs  (LOX, 
TBX,  NLRP1  In  black)  and  SEs  (FS3,  CCL2, 
irCAMf  in  ted)  in  response  to  TNF  a  and  JQ1  po, 
100, 250, 500  nM).  The  mRNA  levels  from  TNF  a  + 
VEH  (10  ng/ml,  3  hr)  treated  ECs  were  set  to 
100%.  Results  displayed  as  the  percent  reduction 
from  maximum.  Error  bars  represent  SEM. 
Representative  results  from  two  Independent  ex 
perlments  are  shown. 

(E)  Gene  tracks  from  Chem  seq  (JQ1)  and  ChIP 
seq  PRD4,  H3K27ac)  data  sets  of  the  SELE  SE 
locus  (rpm/bp)  for  JQ1,  BRD4,  and  H3K27acftDm 
TNF  a^  )  or  TNF  a(+)  stimulated  ECs. 

(F  and  G)  Scatter  plot  of  JQ1  getxxne  widebindino 
levels  on  the  X  ax  Is  compared  to  the  log:  change  in 
BRD4  (F)  or  H3K27ac  (0)  ChIP  seq  signal  on  the  y 
axis.  The  change  In  BRD4  and  l-l3K27ac  signal 
was  determined  comparing  TNF  a  +  JQ1  versus 
TNF  a. 

See  also  Figures  S5  and  S6. 


histone  mark  defining  active  enhancers  that  can  be  used  to  iden¬ 
tify  SEs  in  the  absence  of  specific  enhancer  coactivator  data 
(Hnisz  et  ai.,  201 3;  Whyte  et  ai.,  201 3). 

Large  gains  of  histone  acetylation  were  identified  in  macro¬ 
phages  foilowing  LPS  stimulatton  at  the  Irak2  proinflammatory 
gene  locus  (Figure  S6A).  LPS  treatment  provoked  an  immediate 
increase  in  promoter  and  intergenic  acetylation,  accompanied 
by  an  increase  in  RNA  Pol  II  enrichment  throughout  the  gene 
body.  Globally,  SE  analysis  identified  122  gained  SEs  foilowing 
LPS  treatment  near  proinflammatory  genes  involved  in  cytokine 
signaiing,  ceil  adhesion,  and  chemotaxis  (Rgures  S6B  and  S6C). 
Unlike  inflammatory  ECs,  the  target  genes  induced  and  associ¬ 


ated  with  proinflammatory  SEs  in  macro¬ 
phages  are  iargeiy  distinct  from  those 
found  in  ECs  (Figure  S6D).  Transcription 
at  SE-associated  genes  was  more 
strongly  induced  by  LPS  compared  to 
genes  controlied  by  typical  enhancers, 
as  revealed  by  composite  analysis  of  RNA  Poi  II  enrichment 
and  elongation  (Rgures  S6E  and  S6F).  As  in  activated  ECs, 
BET  bromodomain  inhibition  in  macrophages  preferentially  sup¬ 
pressed  transcription  of  genes  driven  by  proinflammatory  SEs  as 
compared  to  genes  controlled  by  typical  enhancers  (Rgures 
S6G  and  S6H).  The/rak2  locus  provides  an  exemplary  illustration 
of  the  effect  of  BET  bromodomain  inhibition  using  the  structurally 
analogous  l-BET  inhibitor  (Nicodeme  et  al.,  2010),  on  depletion 
of  RNA  Pol  II  enrichment  (Figure  S6A).  Collectively,  these  data 
from  ECs  and  macrophages  demonstrate  that  proinflammatory 
SEs  drive  proinflammatory  gene  activation  in  a  celi  context  spe¬ 
cific  manner. 
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BET  Bromodomain  Inhibition  Suppresses  Leukocyte 
Rolling,  Adhesion,  and  Transmigration  in  Endothelium 
In  response  to  proinflammatory  stimuli,  endotheliaHeukocyte 
interactbns  follow  a  sequential  cascade  involving  ieukocyte 
chenrK3attraction,  their  slow  rolling,  and  subsequent  firm  adhe¬ 
sion  to  ECs,  cuiminating  in  leukocyte  endothelial  transmigration 
into  tissue  (Ley  et  al.,  2007).  To  explore  a  phenotypic  effect  of 
SE  disruption  by  BET  bromodomain  inhibition,  we  tested  the 
functional  effects  of  JQ1  on  leukocyte  rolling  across  TNF-a-acti- 
vated  endothelium  in  vivo.  C57BI/6  mice  were  pretreated  with 
JQ1  (50  mg/kg)  12  hours  before  TNF-a  injection.  As  depicted 
with  intravital  microscopy  of  leukocyte  rolling  in  the  cremaster 
postcapillary  venule,  BET  bromodomain  inhibition  significantly 
reduced  the  leukocyte  rolling  flux  (1 5.8  versus  7.5,  p  <  0.01 ;  Fig¬ 
ure  6A)  and  the  number  of  leukocyte  rdlers/minute  (42.8  versus 
25.14,  p  <  0.05;  data  not  shown),  without  changing  systemic 
white  blood  cell  count  or  shear  stress  (Figure  S7A).  BET  bromo¬ 
domain  inhibition  also  shifted  the  distribution  of  leukocyte  veloc¬ 
ity  and  increased  mean  velocity  (2.89  pm/s  versus  3.91  pm/s,  p  < 
0.01 ),  consistent  with  an  effect  on  E-selectin-mediated  slow  rolF 
ing  (Figures  6B  and  6C). 

Next,  we  examined  firm  adhesion  of  the  human  monocytic  cell 
line  (THP1)  to  activated  ECs  in  vitro.  TNF-a-stimulated  ECs  had 
significantly  increased  numbers  of  attached  THPI  cells  (Fig¬ 
ure  6D).  JQ1  pretreatment  of  ECs  suppressed  THPI  adhesion 
to  TNF-a-activated  ECs  by  70%  (Figures  6D  and  6E).  Similarly, 
siRNA  inhibition  of  BRD4  expression  in  ECs  recapitulated  JQ1  ’s 
effects  on  THPI  adhesion  to  ECs  (Figures  6F  and  6G).  Lastly, 
we  tested  BET  bronrK)domain  inhibitor  effects  on  leukocyte  trans¬ 
migration  in  a  p>arallel-plate  flow  chamber.  TNF-a  stimulatbn  of 
ECs  resulted  In  transmigratbn  of  67%  of  human  neutrophils  (Fig¬ 
ure  6H).  Pretreatment  of  ECs  with  JQ1  prior  to  TNF-a  stimulation 
reduced  neutrophil  transmigration  in  a  concentration-deperbent 
manner  (Figure  6H).  In  kinetb  studies  of  transcription  response  of 
TNF-a-stimulated  ECs,  JQ1  demonstrated  prolonged  inhibitory 
effects  on  expressbn  of  SE-associated  genes  {SELE,  VCAM1, 
CXCL8,  CCL2)  over  48  hr  (Rgures  6I-6L).  These  data  establish 
BET  bromodomain  inhibitbn  as  a  furtctbnal  suppressor  of  the 
phenotypic  features  of  EC  proinflammatory  activatbn. 

BET  Bromodomain  Inhibition  Suppresses 
Atherogenesis  in  Hypercholesterolemic  Mice 
Proinflammatory  activation  of  ECs  is  a  seminal,  early  event  in 
atherogenesis,  a  process  driven  by  vascular  inflammation  that 
also  involves  monocytes/macrophages  and  precedes  athero¬ 
sclerosis  (Cybulsky  et  al.,  2001).  The  transcript bna I  and  func¬ 
tional  effects  of  BET  bromodomain  inhibition  in  ECs  prompted 
us  to  examine  their  role  in  atherogenesis  using  the  well-estab¬ 
lished  low-density  lipoprotein  (LDL)  receptor-deficient  (Ldlr"^~) 
mouse  model.  Vehicle-treated  mice  fed  a  cholesterol-enriched 
diet  (1 0  weeks)  developed  atheroscbrosis,  as  measured  by  oil 
red  O  staining  (Figure  7A).  Orxie-daily  JQ1  treatment  (50  mg/ 
kg)  reduced  aortic  plaque  area  by  40%  (Figure  7A).  Notably, 
there  was  no  difference  in  LDL  artd  high-density  lipoprotein 
cholesterol  between  the  vehicle  and  JQ1  treatment  groups 
(168  versus  164  mg/dL  and  56  versus  59  mg/dL,  respectively, 
total  cholesterol  911  versus  1,349  mg/dl;  see  also  Table  S6.). 
Early  atherosclerotic  lesions  are  comprised  of  macrophages 


(98%-99%)  with  tesser  amounts  of  T  lymphocytes  (1%-2%). 
Mac-3  staining  demonstrated  that  JQ1  treatment  significantly 
lowered  total  macrophage  staining  area,  CD4-positive  T  lympho¬ 
cytes,  and  levels  of  VCAM1  protein  (Figures  7B-7D).  Oil  red  O 
staining  of  en  face  thoracoabdominal  aortas  revealed  decreased 
atherosclerotic  plaque  beyond  the  aortic  root  (Figures  7E  arb 
7F).  Soluble  VCAM1  and  ICAM1  bvels  were  also  significantly 
reduced  in  JQ1 -treated  animals  compared  to  vehicle  (Figures 
S7B  arb  S7C),  suggesting  an  effect  of  BET  bromodomain  inhibi¬ 
tion  on  systemb  proinflammatory  activation  in  vivo.  We  next 
tested  whether  BET  bronxDdomain  inhibitbn  mitigated  the  acti¬ 
vation  of  proad  hesbn  prathways  in  aortic  erbothelium,  which  oc¬ 
curs  during  the  first  1 0  weeks  of  exposure  to  an  atherogenb  diet. 
In  ex  vivo  aortic  adhesbn  assays,  the  aortas  harvested  from  an¬ 
imals  (6  \Meekson  dbt)  treated  with  JQ1  supported  less  adhesion 
of  fluorescently  labeled  monocytes  (Rgures  S7D  and  S7E). 
These  data  demonstrate  BET  inhibition  significantly  decreased 
atherogenesis  and  accumulation  of  inflammatory  cells  in  a 
welFcharacterized  murine  model  of  atherosclerosis. 

DISCUSSION 

The  inflammatory  response  underlies  numerous  chronb  dis¬ 
eases.  NF-kB  is  a  master  regulatory  transcriptbn  factor  in 
several  domirtant  inflammatory  signaling  cascades,  which  coop¬ 
erate  with  chromatin-associated  regulatory  comptexes  to  direct 
inflammatory  transcription  (Natoli,  2009).  Enhancer-bound  NF- 
kB  arises  following  nuclear  translocation  in  a  manner  influenced 
by  pioneer  transcription  factors  (Kaikkonen  et  al.,  2013;  Natoli, 
2009;  Ostuni  et  al.,  2013)  and  a  preestablished  topwlogy  that 
places  distal  enhancer  regions  and  target  genes  in  spratial  prox¬ 
imity  (Jin  et  al.,  2013).  Here,  we  explore  the  role  of  chromatin  in 
terminal  signal  transduction  from  NF-kB  to  RNA  polymerase, 
specifically  at  massive  regulatory  regions. 

Super  or  stretch  enhancers  represent  less  than  5%  of  the  en¬ 
hancers  in  a  cell,  yet  they  use  almost  half  of  all  enhancer  coac¬ 
tivator  proteins  (Loven  et  al.,  201 3;  Parker  et  al.,  201 3;  Whyte 
et  al.,  2013)  arb  are  highly  transcribed,  producing  large  amounts 
of  enhancer  RNA  that  may  itself  facilitate  target  gene  activation 
(Hnisz  et  al.,  2013;  Kaikkonen  etal.,  2013).  The  SE  landscape  is 
remarkably  cell-type  specific,  driving  expression  of  the  genes 
that  define  and  maintain  cell  identity  in  different  tissues  arb 
cell  lineages.  The  present  study  demonstrates  that  NF-kB  en¬ 
gages  most  endothelial  enhancers  following  proinflammatory 
activation,  yet  significant  BRD4  recmitment  to  form  de  novo 
SEs  is  restricted  to  a  subset  of  these  enhancer  regions.  NF- 
KB-directed  SE  formation  causes  gbbal  reorganizatbn  of  the 
BRD4  SE  larbscape  arb  induces  the  transcriptbn  of  many  ca- 
rxinical  proinflammatory  endothelial  genes.  Whereas  previous 
studies  have  studied  the  importance  of  SEs  in  the  maintenance 
of  cell  identity  (Whyte  et  al.,  201 3),  here  we  describe  de  novo  SE 
formation  as  a  mechanism  by  whbh  stimulus-coupled  master 
regulatory  transcriptbn  factors  such  as  NF-kB  can  coordinate 
a  rapid  transcriptbnal  response  that  drives  a  dynamic  change 
in  cell  state. 

This  study  of  kinetb  transcriptional  response  during  the  in¬ 
flammatory  cell  state  transition  in  ECs  unexpectedly  bentified 
a  rapb  loss  of  SEs  upon  cytokine  stimulatbn.  As  orchestrated 
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Rgure  6.  Phenotypic  Consequences  of  BET 
Bromodomain  Inhibition  in  Endothelium 

(A)  Intravltal  microscopy  Image  (left)  and  bar  plot 
quantification  (right)  of  leukocyte  flux  fraction  In  the 
cremaster  postcaplllary  venule  after  TNF  <x  (2  hr, 
n  =  7/group)  in  VEH  orJQI  treated  samples.  En'or 
bars  nepresent  SEM.  The  statistical  significance 
of  the  difference  between  JQ1  and  VEH  treated 
samples  was  determined  using  a  two  tailed  t  test. 
*p  <  0.05. 

(B)  Velocity  distribution  of  leukocytes  measured  In 
(A). 


01 


K 


^)  Bar  plot  showing  mean  leukocyte  velocity  in 
cremaster  postcaplllary  venule  In  TNF  a{+) 
animals  ±  BET  bromodomain  liTilbItlon.  Error  bars 
represent  SEM. 

P  and  F)  Representative  fiuoresoenoe  microscopy 
images  showing  adhesion  of  calcein  labeled  THP1 
cells  to  (D)  ECs  pretreated  with  JQ1  then  activated 
with  TNFa  (4  hr)  as  well  as  (F)  TNF  a  treated  ECs 
after  sIRNA  knockdown  of  BRDA. 

(E  and  Q)  Bar  plots  showing  quantification  of 
fluorescence  from  (D)  and  (F). 

(H)  Bar  plots  showing  quantification  of  trans 
migrating  neutrophils  on  TNF  a  activated  EC 
monolayers.  Results  pooled  from  three  Indepen 
dent  experiments.  Data  represent  mean  ±  SEM. 
The  statistical  significance  of  the  difference  be 
tween  JQ1  and  VEH  treated  samples  was  deter 
mined  using  a  two  tailed  t  test.  *p  <  0.05. 

(1  L)  Line  plots  of  mRNA  levels  (qRT  PCR)  for  (I) 
SELE,  (J)  VCAM7,  (K)  CXCL8,  and  (L)  CXX2 
measured  after  stimulation  of  ECs  with  TNF  o 
(1 2.5  ng/ml;  1 , 3. 8.  24,  48  hr)  ±  JQ1  (500  nM). 

The  statistical  significance  of  the  difference  be 
tween  samples  was  determined  using  a  two  tailed 
t  test.  *p  <  0.05  in  TNF  o(+)  versus  TNF  o(  ); 
ftp  <  0.05  in  JQ1  versus  VEH.  Data  represent 
mean  ±  SEM  of  fold  change  versus  0  hr.  See  also 
Figure  S7. 
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by  the  master  regulatory  inflammatory  transcriptbn  factor 
NF-kB,  SE  formation  comes  at  the  immediate  expense  of  SEs 
associated  with  active  transcription  of  genes  in  unstimulated 


for  enhancer  output, 
thereby  driving  rapid 
transcription. 


ECs,  including  targets  involved  in  cell 
specification  and  noninflammatory  cell 
states.  Many  of  these  tost  SE  associated 
genes,  such  as  SOX18,  will  be  down- 
regulated.  Curiously,  we  find  persistent 
marks  of  open,  active  euchromatin  at 
the  SOX18  locus  {H3K27ac),  yet  BRD4 
depletion  is  associated  with  toss  of  RNA 
pol  II  transcription.  In  studies  of  dynamic 
cell  state  transitions,  as  here,  we  find 
that  experimental  measurements  of  ca¬ 
nonical  enhancer  marks  (i.e.,  H3K27ac, 
H3K4me1/2)  are  inferior  to  assessments 
of  enhancer  coactivators  (i.e.,  BRD4)  for 
characterizing  rapid  changes  in  enhancer 
structure  and  function.  As  such,  BRD4 
may  be  viewed  and  used  as  a  rheostat 
converting  typical  enhancers  into  SEs, 
and  robust  induction  of  inflammatory 
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In  our  previous  studies,  genes  with  the  highest  occupancy 
of  BRD4  at  their  proximal  SEs  were  the  most  selectively 
downregulated  by  BET  bromodomain  inhibitbn  (Chapuy 
et  al.,  2013;  Loven  et  al.,  2013).  However,  these  studies 
characterized  effects  on  cells  at  ground  state,  where  stable, 
preestablished  SEs  predominate.  During  cell  state  transi¬ 
tions,  such  as  the  present  study  of  EC  and  macrophage 
activation,  we  observe  potent  and  selective  effects  on  up- 
regulated  genes  associated  with  de  novo  SEs.  Preferential 
disruption  of  dynamic  SEs  by  BET  bromodomain  inhibition 
abrogates  the  induction  of  inflammatory  transcription.  The 
direct  relationship  between  BRD4  enrichment  and  transcrip¬ 
tion  suggests  that  modulating  BRD4  levels  at  enhancers 
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Figure  7.  BET  Bromodomain  Inhibition  Sup 
presses  Atherogenesis  in  Ldlr~'~  Mice 
(A  D)  Photomicrographs  of  aortic  root  sections 
from  Ldir  '  animals  treated  with  VEH  or  JQ1 
stained  tor  (A)  oil  rad  O,  (B)  Mac  3,  (C)  CD  4,  or  (D) 
VCAM1 .  Quantification  of  staining  is  shown  below. 
Results  represent  mean  ±  SEM.  The  statistical 
significance  of  the  difference  between  JQ1  and 
VEH  treated  samples  was  determined  usbig  a 
two  tailed  t  test.  *p  =  0.002  for  (A);  atxl  *p  <  0.05  for 
(B)  (D). 

(E)  Oil  red  O  staining  of  en  face  aortas  prepared 
from  cohort  In  (A)  (D). 

(F)  Quantification  of  lesion  area  (%)  between  VEH 
and  JQ1  treated  en  face  aortas.  The  statistical 
significance  of  the  difference  between  JQ1  and 
VEH  treated  samples  was  determined  using  a 
two  tailed  t  test.  *p  <  0.05. 

See  also  Figure  S7. 
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may  influence  the  pathogenesis  of 
inflammatory  diseases.  Changes  in 
eRNA  levels  by  NF-KB-directed  BRD4 
SE  formation  may  also  be  an  impiortant 
factor  in  proinflammatory  transcription 
in  ECs,  but  await  future  study.  Further 
granuiarity  on  dynamic  enhancer  re- 
modeiing  wiii  accompany  progress  in 
geixime-wide  enhancer  detectbn  and 
assignment. 

In  many  disease  settings,  the  degree 
of  the  host  infiammatory  response  is  a 
key  determinant  of  severity  (Medzhitov 
et  al.,  2012).  Here  we  show  small 
molecule  BET  bromodomain  inhibition 
(JQ1)  significantiy  attenuated  endotheiiai 
activation  during  acute  infiammation 
in  vitro  and  ex  vivo,  as  reveaied  by  sup- 
pressbn  of  TNF-a  induced  ieukocyte 
rdiing,  adhesion  and  transmigration. 
Finaily,  in  a  hyperchoiesteroiemic  murine 
nrodei  of  atheroscierosis,  in  whbh  EC 
proinfiammatory  activation  is  a  seminai 
eariy  event,  10-week  treatment  with 
JQ1  suppressed  cardinai  histopathobgic 
features  of  atherogenesis.  These  data  in  vascular  endothelial 
activation  estabiish  a  critbal  and  eariy  roie  for  BET  bromodo- 
mains  in  dynamb  enhancer  renrKXleling.  They  describe  a 
mechanism  for  rapid  infiammatory  gene  activation  by  NF-kB- 
mediated  formation  of  SEs.  Taken  in  the  context  of  prior 
research  in  macrophage  activation,  spermatogenesis,  and 
myocyte  hypertrophy  (Anand  et  al.,  2013;  Delmore  et  al., 
2011;  Matzuk  et  al.,  2012;  Nbodeme  et  al.,  2010),  these  data 
support  a  model  where  localizatbn  of  BET  bromodomains 
to  chromatin  facilitates  cell  state  transitbns.  The  existerx^e  of 
BET  bromodomain  inhibitors  provides,  then,  a  broad  oppor¬ 
tunity  for  infiammatory  gene  control  through  modulatbn  of 
chromatin  structure  and  functbn. 
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EXPERIMENTAL  PROCEDURES 
Animal  Models 

LDL  receptor  knockout  mice  (4  weeks  old)  on  a  C57BI/<6  background  were 
purchased  from  Jackson  Laboratories.  Mice  were  fed  an  atherogenic  diet 
(Clinton/Cybulsky  Rodent  Diet,  D12108  with  1.25%  cholesterol,  Research 
Diets)  for  10  weeks.  While  on  the  diet,  the  animals  were  treated  with  vehicle 
(DMSO)  or  JQI  (SO  mg/kg)  by  intraperltoneal  Injection,  once  dally  (N  =  10/ 
group).  Oil  red  O  staining  was  used  to  quantify  atherosclerotic  plaque  lesion 
area  In  the  aortic  root.  Macrophage  (Mac  3)  and  T  lymphocyte  (CD  4)accumu 
latlon  was  assessed  in  the  aortic  root,  and  the  total  staining  area  was  analyzed 
using  computer  assisted  Imaging  analysis.  For  CD  4  cells,  total  cell  numbers 
were  counted.  Whole  aorta  from  the  left  subclavian  artery  to  the  iliacbifurcation, 
was  used  for  en  face  preparation  and  stained  with  oil  red  O.  All  protocols  con 
ceming  animal  use  were  approved  by  the  Harvard  Medical  School  Institutional 
/Vnimal  Care  and  Use  Committee  and  conducted  In  accordance  with  the  NIH 
Guide  forthe  Care  and  Use  of  Laboratory  Animals.  All  studies  were  performed 
In  C57BI/6J  mice  (Jackson  Laboratories),  maintained  in  a  pathogen  free  facility 
with  standard  light/dark  cycling  and  ad  libitum  access  to  food  and  water. 

Reagents  and  Cell  Culture 

ECs  from  pooled  human  umbilical  cords  were  cultured  In  M199  medium 
supplemented  with  20%  fetal  bovine  serum  (FBS),  0.1%  heparin,  50  pgAnI 
endothelial  cell  growth  factor  (Biomedical  Technologies),  peniclllin/strepto 
mycin  on  gelatin  coated  tissue  culture  plates.  U937  cells  (ATCC)  were  main 
talned  In  RPMI  with  10%  FBS  and  antibiotics.  For  transendothellal  migration 
(see  below),  human  umbilical  vein  ECs  (subculture  2)  were  grown  on  fibro 
nectin  coated  glass  coverslips  (5  mg/ml;  BD  Biosciences)  and  treated  with 
JQI  (500  nM)  or  vehicle  (DMSO)  for  1  hr  before  TNF  a  stimulation  (10  ng/ml, 
4  hr).  Recombinant  human  TNF  a  was  obtained  from  PeproTech.  JQI  was  dis 
solved  In  DMSO  at  a  concentration  of  50  mg/ml.  Working  stocks  of  JQI  were 
prepared  by  diluting  1 :10  In  10%  beta  cyclodextrin  solution  (Filippakopoulos 
et  al.,  2010).  /Vnimals  were  treated  at  50  mg/kg  once  dally  by  Intraperltoneal 
Injection. 

ChIP-Seq  and  Data  Analysis 

Chip  was  performed  in  ECs  in  the  presence  or  absence  of  TNF  a  (1  hr,  25  ng/ 
ml)  and  JQI  (3  hr  pretreatment,  500  nM).  Specific  antibodies  ard  detailed 
methods  are  described  In  the  Supplemental  Experimental  Procedures.  ChIP 
was  carried  out  as  described  elsewhere  (Lov^  et  al.,  2013).  All  ChIP  seq 
data  sets  were  aligned  using  Bowtie  (version  1.0.0;  Langmead  et  al.,  2009) 
to  build  version  NCBI36/HG18  of  the  human  gerxxne  or  build  version 
NCBI37/MM9  of  the  mouse  genome.  Enhancers  and  super  enharcers  were 
mapped  as  described  elsewhere  (Lov^  et  al.,  2013).  /Additional  details  are 
provided  in  the  Supplemental  Experimental  Procedures. 

ACCESSION  NUMBERS 

The  GEO  fittp://www.ncblnlm.nih.gov/geo)  accession  numbers  for  all  ChIP 
seq  and  Chem  seq  data  (bxiudkig  microarray  dateO  and  aligned  and  raw 
data  reported  In  this  work  are  GSE53998  and  GSE54000,  respectively. 

SUPPLEMENTAL  INFORMATION 

Supplemental  Information  includes  Supplemental  Experimental  Procedures, 
seven  figures,  and  six  tables  artd  can  be  found  with  this  article  online  at 
http;//dx.doi.org/1 0. 101  €/j.molcel.2014.08.024. 
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Summary  (148  words): 


Medulloblastoma  is  a  highly  malignant  paediatric  brain  tumour,  often  inflicting 
devastating  consequences  on  the  developing  child.  Genomic  studies  have  revealed  four 
transcriptionally  distinct  molecular  subgroups  with  divergent  biology  and  clinical 
behaviour.  An  understanding  of  the  regulatory  circuitry  governing  the  transcriptional 
landscapes  of  medulloblastoma  subgroups,  and  how  this  relates  to  their  respective 
developmental  origins,  is  currently  lacking.  Using  H3K27ac  and  BRD4  ChIP-Seq, 
coupled  with  tissue-matched  DNA  methylation  and  transcriptome  data,  we  describe  the 
active  c/s-regulatory  landscape  across  28  primary  medulloblastoma  specimens.  Analysis 
of  differentially  regulated  enhancers  and  super-enhancers  reinforced  inter-subgroup 
heterogeneity  and  revealed  clinically  relevant  insights  into  oncogenic  TGPp  signaling  in 
Group  3.  Computational  reconstruction  of  core  regulatory  circuitry  identified  a  master  set 
of  transcription  factors  responsible  for  subgroup  divergence  and  revealed  candidate 
cellular  origins  for  Group  4.  The  integrated  analysis  of  c/s-regulatory  elements  in  primary 
human  tumour  samples  reveals  insights  into  c/s-regulatory  architecture,  unrecognized 
dependencies,  and  cellular  origins. 


Introduction  (465  words): 


Medulloblastoma  is  a  highly  malignant  paediatric  brain  tumour  classified  into  four 
biologically  and  clinically  distinct  molecular  subgroups^'^.  Transcriptional  diversity 
underlying  WNT,  SHH,  Group  3,  and  Group  4  subgroup  medulloblastoma  is  partially 
explained  by  active  and  discriminatory  signaling  pathways,  such  as  the  Wingless/WNT 
and  Sonic  hedgehog/SHH  developmental  cascades  inherent  to  WNT  and  SHH 
medulloblastomas,  respectively.  Recurrent,  somatically  altered  driver  genes  such  as 
MYC  (Group  3),  KDM6A  (Group  4),  the  recently  implicated  GFI1/GFI1B  (Group  3  and 
Group  4),  and  others  contribute  to  further  diversity  between  medulloblastoma 
subgroups'*'®.  Whereas  distinct  cellular  origins  have  been  experimentally  substantiated 
for  WNT^  and  SHH®'*°  tumours,  the  origins  of  Group  3  and  Group  4  medulloblastoma 
remain  unknown.  Underscoring  the  need  for  new  subgroup-specific  therapeutic  insights, 
the  present  clinical  approach  to  medulloblastoma  involves  invasive  surgery,  cranio¬ 
spinal  radiation  and  cytotoxicity,  together  associated  with  profound  morbidity  in  the 
developing  child. 

Recent  next-generation  sequencing  studies  of  medulloblastoma  have  defined  recurrently 
mutated  genes  and  pathways,  the  proportion  of  cases  affected  by  such  alterations,  and 
their  respective  subgroup  distribution'*  **'*'*.  The  bulk  of  these  efforts  have  thus  far  been 
focused  on  somatic,  DNA-level  genomic  alterations,  especially  non-synonymous  single 
nucleotide  variants,  indels,  and  focal  copy-number  aberrations®.  Recurrent  targeting  of 
genes  involved  in  chromatin  modification  has  been  the  most  consistent  theme  to  emerge 
from  these  studies**  *®  *®,  strongly  suggesting  deregulation  of  the  epigenome  as  a  critical 
step  during  medulloblastoma  pathogenesis.  However,  this  hypothesis  has  yet  to  be 
substantiated  and  knowledge  pertaining  to  how  the  medulloblastoma  epigenome 
influences  subgroup-specific  transcriptional  programs  remains  in  its  infancy.  Recent 
analysis  of  DNA  cytidine  methylation  in  medulloblastoma  corroborated  transcriptional 
differences  in  subgroups  and  enumerated  putative  master  regulatory  factors*^.  Still,  a 
detailed  analysis  of  the  c/s-regulatory  epigenome  of  this  disease  has  not  been 
undertaken. 

Enhancers  are  c/s-acting  regulatory  elements  that  serve  as  sites  of  recruitment  for 
transcription  factors  and  chromatin-associated  regulatory  complexes,  which  together 
signal  to  RNA  polymerase  in  control  of  target  gene  expression*®.  Massive  catalogues  of 
genome-wide  enhancers  have  been  inferred  and  published  by  large  consortia  such  as 
ENCODE*®'^®  and  Roadmap^*,  dramatically  advancing  our  understanding  of  enhancer- 
gene  regulation  across  a  comprehensive  spectrum  of  cell  lines  and  tissues  from  different 
species.  These  resources  empower  our  understanding  of  the  complex  cartography  of  the 
human  regulatory  landscape,  provide  testable  hypotheses  regarding  disease-risk 
association,  contribute  evolutionary  inferences,  and  establish  facile  analytical 
techniques.  To  deeply  characterize  the  active  c/s-regulatory  circuitry  of  a  single  disease 
entity,  here  medulloblastoma,  we  performed  high-resolution  chromatin 
immunoprecipitation  with  sequencing  (ChIP-Seq)  for  active  enhancers  (H3K27ac)  in  28 
primary  tumour  specimens  and  three  established  cell  lines,  collectively  accounting  for 
known  molecular  subgroups.  Our  approach  to  studying  enhancers  genome-wide  in  a 
large  set  of  primary  tissue  samples  led  to  a  regulatory  explanation  for  subgroup 
transcriptional  diversity,  previously  unrecognized  subgroup-specific  dependencies  and 
firm  insights  into  medulloblastoma  cellular  origins,  in  particular  for  the  poorly 
characterized  and  aggressive  Group  3  and  Group  4  subgroups. 


Results:  (2,491  words): 

The  enhancer  landscape  of  primary  medulloblastoma 

Recent  publications  of  large-scale  efforts  to  annotate  active  regulatory  elements 
genome-wide  in  human  tissues  {e.g.  through  DNase  I  hypersensitivity,  H3K27ac  and 
BRD4  ChlP-seq)^®’^\  have  focused  on  immortalized  or  malignant  cell  lines  and  normal 
human  tissues  for  cataloguing  active  enhancers.  Discrete  disease  entities  have  been 
under-represented  in  these  comprehensive  surveys,  such  as  medulloblastoma,  which 
has  been  the  subject  of  only  a  single  long-term  culture  cell  line  (D721;  first  reported  in 
1997)  included  amongst  125  cell  types  initially  studied  by  ENCODE^”.  Further,  cancer 
cell  lines  often  exhibit  drastic  genomic  and  transcriptional  divergence  from  their 
corresponding  primary  tumour  tissues.  This  is  exemplified  in  Non-Hodgkin’s  lymphoma 
where  our  prior  epigenomic  analyses  identified  greater  likeness  between  primary  tumour 
samples  and  normal  lymphoid  tissues  than  between  tumours  and  cell  lines^^.  Given  the 
apparent  limitations  of  using  cell  lines  to  faithfully  study  the  tumour  epigenome,  and  the 
recognized  subgroup-dependent  heterogeneity  of  medulloblastoma,  we  collected  a 
series  of  28  treatment-naive,  fresh-frozen  medulloblastoma  specimens  for  studying  the 
active  enhancer  landscape  by  H3K27ac  ChIP-Seq  (Figure  1a,b,  Extended  Data  Figure 
1a,b). 

The  cohort  was  selected  to  be  inclusive  of  all  four  medulloblastoma  subgroups 
(Supplemental  Table  SI;  WNT,  n=3,  SHH,  n=5,  Group  3,  n=9,  Group  4,  n=11).  Three 
additional  Group  3  cell  lines  (MED8A,  D425,  and  HD-MB03)  were  also  included  in  our 
experimental  workflow.  Using  MACS^^  to  identify  significantly  enriched  H3K27ac  peaks, 
we  inferred  78,516  medulloblastoma  enhancers,  which  mainly  (~80%)  covered  introns 
and  intergenic  regions  (Extended  Data  Figure  1b).  Parallel  ChIP-Seq  was  performed  for 
Bromodomain  Containing  4  (BRD4),  an  enhancer-associated  transcriptional 

coactivator^^’^'*,  in  27/31  cases  (Figure  1c).  Enrichment  of  H3K27ac  and  BRD4  ChIP-Seq 
signals  was  highly  correlated  at  putative  enhancer  loci  (Pearson  correlation,  r=0.949), 
confirming  that  the  selected  regions  are  indeed  active  enhancers  (Figure  Ic)^^’^'*.  In 
contrast,  regions  enriched  for  H3K27ac  were  strongly  anti-correlated  with  DNA 
methylation  (Pearson  correlation,  r=-0.577;  Figure  Id).  Finally,  analyzing  strand-specific 
RNA-Seq  data  generated  from  the  same  medulloblastomas  subjected  to  ChIP-Seq,  we 
observed  short,  unspliced,  bidirectional  RNA  transcripts  overlapping  active  enhancers 
(Figure  1e),  in  accordance  with  recently  described  enhancer  RNAs  (eRNAs),  known  also 
to  arise  from  active  enhancers^^.  Comparison  of  predicted  medulloblastoma  enhancers 
with  those  reported  using  analogous  methods  employed  by  the  ENCODE  and  Roadmap 
Epigenomics  Projects  revealed  19,850  novel  regulatory  regions,  indicative  of  potentially 
cerebellar  cell  type-  or  medulloblastoma-specific  enhancers  in  our  dataset  (Figure  If,  g). 
Importantly,  primary  medulloblastoma  enhancer  landscapes  exhibited  poor  overlap  and 
correlation  with  those  generated  from  medulloblastoma  cell  lines  (Extended  Data  Figure 
1c,  d),  further  emphasizing  the  importance  of  studying  the  epigenome  in  primary 
tumours. 

Since  medulloblastoma  subgroups  were  first  described  based  on  their  transcriptional 
diversity^®,  we  sought  to  explore  subgroup-specific  enhancers  to  gain  insight  into  the  cis- 
regulatory  landscape  defining  subgroup  identity.  ANOVA  was  used  to  identify  sets  of 
enhancers  differing  according  to  known  molecular  subgroup,  revealing  20,406 
differentially  regulated  enhancers  (26%  of  all  inferred  medulloblastoma  enhancers; 
Figure  2a,  b).  The  remaining  74%  (n=58,110)  displayed  activity  across  all  groups 


suggesting  a  general  role  in  medulloblastoma  or  cerebellar  identity  (Figure  2a; 
Supplemental  Table  S2).  K-means  clustering  of  differentially  regulated  enhancers 
delineated  six  distinct  medulloblastoma  enhancer  classes,  including  one  for  each 
subgroup  (i.e.  WNT,  SHH,  Group  3  and  Group  4)  as  well  as  WNT-SHH  and  Group  3- 
Group  4  shared  classes  (Figure  2b,  c).  Group  3  and  Group  4  subgroups  are  known  to 
exhibit  some  degree  of  transcriptional  similarity^^'^^,  consistent  with  the  enhancer 
clustering  results.  In  contrast,  WNT  and  SHH  subgroups  tend  to  be  mostly  dissimilar 
from  a  transcriptional  perspective,  and  thus  a  common  subset  of  shared  enhancers 
between  these  groups  was  unexpected  and  intriguing. 

We  next  sought  to  assign  enhancer  elements  to  target  genes,  a  process  typically 
hindered  by  the  fact  that  the  majority  of  enhancer/promoter  interactions  occur  over 
extensive  genomic  distances^”.  To  overcome  challenges  in  enhancer/gene  assignment, 
we  leveraged  sample-matched  RNA-Seq  gene  expression  data  to  identify  putative 
enhancer/gene  interactions  that  are  (i)  contained  in  the  same  topologically  associated 
domain  (TAD^^)  and  (ii)  exhibit  strong  positive  correlations  between  enhancer  H3K27ac 
levels  and  gene  expression  (p>  0.6  and  FDR  <0.05,  Extended  Data  Figures  2  &  3). 
TADs  are  megabase  scale  genomic  regions  of  interacting  chromatin  that  form  a 
fundamental  unit  of  genome  structure.  This  approach  assigned  8,775  enhancers 
(including  43%  of  all  differential  enhancers)  to  at  least  one  protein-coding  target  gene 
(Supplemental  Table  S3).  The  majority  (44%)  of  inferred  target  genes  were  assigned  to 
a  single  enhancer,  but  in  many  cases,  several  enhancers  were  predicted  to  converge  on 
the  regulation  of  a  single  gene  (Figure  2d).  Likewise,  73%  of  enhancers  were  assigned 
to  only  a  single  gene  target  and  rarely  was  a  given  enhancer  assigned  to  more  than  two 
candidate  genes  (Figure  2e).  Compelling  subgroup-related  diversity  with  respect  to 
inferred  enhancer-target  gene  regulation  was  prevalent  in  our  dataset  (Figure  2f-i),  with 
numerous  genes  exhibiting  convergent  regulation  by  distinct  subgroup-specific 
differential  enhancer  loci.  For  example,  we  identified  alternative  subgroup-specific 
enhancers  predicted  to  regulate  known  oncogenes,  including  WNT-specific  and  SHH- 
specific  enhancers  inferred  to  target  ALK,  and  WNT-specific  and  Group  3-specific 
enhancers  inferred  to  target  MYC  (Figure  2g-i;  Extended  Data  Figure  2j).  These  data 
provide  a  rational,  computationally  robust  approach  for  assigning  medulloblastoma 
enhancers  to  their  potential  targets  and  underscore  the  apparent  complexity  inherent  to 
enhancer-gene  regulation  across  medulloblastoma  subgroups. 

Subgroup-specific  enhancers  reveai  aberrant  TGFfi  signaiing  in  Group  3 

Group  3  and  Group  4  medulloblastoma  remain  the  least  well  understood  subgroups  of 
the  disease,  despite  collectively  accounting  for  -60%  of  all  diagnoses^’^\  Group  3 
patients  have  a  particularly  grave  prognosis,  suffer  commonly  from  metastatic  disease, 
and  fall  ill  in  infancy  or  early  childhood.  Group  4  disease  affects  children  of  all  ages  and 
prognoses  have  been  demonstrated  to  vary  according  to  metastatic  status  and  the 
presence  of  certain  cytogenetic  anomalies^^.  In  contrast  to  WNT  patients  (who  almost 
universally  survive  current  treatment  regimens),  and  SHH  patients  (who  represent 
rational  candidates  for  SHH  pathway  inhibitors  such  as  SMO  inhibitors),  mechanism- 
based  treatment  options  remain  scarce  for  Group  3  and  Group  4  patients.  As  such, 
additional  insights  are  needed  to  direct  future  mechanistic  and  translational  research. 

We  therefore  applied  differential  enhancer  analysis  to  identify  divergent  regulatory 
programs  between  Group  3  and  Group  4  tumours.  First,  we  validated  that  Group  3  or 
Group  4  differential  enhancer  target  genes  showed  reciprocal  patterns  of  H3K27ac 


enrichment.  We  ranked  the  top  1,000  enhancers  (by  H3K27ac  signal)  in  either  Group  3 
or  Group  4  by  fold-change  in  acetylation  and  found  a  strong  leading  edge  enrichment  of 
Group  3-  and  Group  4-specific  target  genes  (Figure  3a).  Functional  pathway  analysis 
(see  Supplementary  Methods)  performed  on  Group  3  and  Group  4  enhancer-gene  target 
assignments  identified  prominent  neuronal  gene  sets  enriched  in  both  subgroups  (Figure 
3b,  c),  consistent  with  published  transcriptional  studies^^'^^  and  validating  the 
computational  approach  implemented  to  assign  enhancers  to  target  genes.  Neuronal 
development  driven  by  transcriptional  regulators  dominated  the  Group  4  functional 
annotation,  whereas  Group  3  enhancer  target  genes  prominently  included  thematic 
pathways  associated  with  TGFp  signaling  (Figure  3b, c;  Extended  Data  Figure  4a). 
Group  3-specific  enhancers  included  the  TGFp  family  type  I  and  II  membrane  receptors 
{ACVR2A  and  TGFBR1)  as  target  genes  (Figure  3c),  whereas  enhancers  regulating 
SMAD5,  TGFB1,  and  TGFB3  showed  equivalent  acetylation  between  Group  3  and 
Group  4.  Overall,  components  of  the  TGFp  signaling  pathway  showed  a  strong 
enrichment  for  enhancer  regulation  in  Group  3  (Figure  3dX  Notably,  we  uncovered  a 
~450kb  focal  amplification  at  the  ACVR2A  locus  in  one  Group  3  sample  that 
encompassed  both  the  gene  and  the  upstream  enhancer  regions  (Extended  Data  Figure 
4b).  These  data,  combined  with  our  prior  observations  that  TGFp  receptor  genes  are 
recurrently  amplified  at  low  frequency  in  Group  3^,  further  implicate  TGFp  signaling  as  a 
putative  oncogenic  driver  in  Group  3  medulloblastoma. 

Medulloblastoma  super-enhancers  define  subgroup  identity 

The  large.  Group  3-specific  enhancer  clusters  adjacent  to  genes  defining  the  TGFp 
signaling  pathway  prompted  a  consideration  that  super-enhancers  (SEs),  broad  spatially 
co-localized  enhancer  domains^^'^®,  might  play  an  essential  role  in  establishing 
subgroup-specific  identity.  SEs  are  established  by  cell  state-defining  transcription  factors 
(TFs)  and  TFs  at  the  termini  of  signaling  pathways®'*®®.  In  multiple  tumour  types,  SEs 
have  been  shown  to  drive  oncogenes,  genes  required  for  maintenance  of  tumour  cell 
identity,  and  genes  associated  with  cell  type-specific  functions.  We  therefore  next 
undertook  a  systematic  mapping  of  SEs  across  all  28  medulloblastoma  samples  (Figure 
4a).  Massive  (>50kb)  SE  domains  were  identified  at  the  cerebellar-specific  TFs,  ZIC1 
and  (Figure  4b,  c),  and  at  70%  of  a  queried  set  of  established  medulloblastoma 

driver/signature  genes,  including  GLI2,  MYC,  OTX2  and  others®  (Extended  Data  Figure 
5a). 

To  identify  subgroup-specific  SEs,  we  ranked  average  H3K27ac  enrichment  across  all 
samples  from  a  given  subgroup  at  the  union  of  all  enhancer  regions  identified  in  that 
subgroup.  Subgroup  SEs  were  identified  from  this  meta  H3K27ac  signal  using  previously 
established  methods®®,  resulting  in  -3,000  distinct  SE  containing  loci  with  -600-1,100 
SEs  identified  per  subgroup  (Figure  4a,  Supplemental  Table  S4).  Compared  to  typical 
enhancers,  subgroup  SEs  showed  higher  occupancy  of  BRD4  and  greater  enhancer 
signal  dynamic  range  between  subgroups  (Extended  Data  Figure  5b-d).  Targets  of 
differential  enhancers  contained  within  SEs  (i.e.  SE  target  genes)  included  a  large 
fraction  of  established  medulloblastoma  signature  genes  (32%;  Supplemental  Table  S3), 
as  well  as  novel  candidates,  including  NKD1/NKD2  (WNT  subgroup),  PONT  (SHH 
subgroup),  HLX  (Group  3),  and  SNCAIP  (Group  4)  (Figure  4d-f).  Medulloblastoma  SEs 
were  inferred  to  regulate  known  Cancer  Gene  Census  genes,  including  ALK  in  WNT, 
SMO  and  NTRK3  in  SHH,  LM01,  LM02,  and  MYC  in  Group  3,  and  ETV4  and  PAX5  in 
Group  4,  among  others  (Supplemental  Table  S3).  Furthermore,  several  actionable,  SE- 


regulated  genes  were  revealed  in  our  analysis  including  kinases  (NTRK1,  SGK1)  and 
chromatin  modifying  enzymes  (PNMT,  HDAC4)  (Supplemental  Table  S5). 


Rank  transformation  of  subgroup  SE  loci  across  all  samples  enabled  a  systematic 
identification  of  SEs  displaying  either  conserved  SE  activity  across  samples,  or  highly 
subgroup-specific  patterns  of  activity  (Figure  4b-f,  Supplemental  Table  S6).  Unbiased 
hierarchical  clustering  of  SEs  across  all  primary  medulloblastoma  samples  was  sufficient 
to  recapitulate  transcriptional  subgroupings  using  no  prior  knowledge  of  subgroup  status, 
suggesting  that  SEs  might  play  a  role  in  driving  subgroup-specific  identity  (Figure  4a).  As 
shown  with  all  enhancer  elements  (Extended  Data  Figure  1c,d),  SEs  from  established 
Group  3  medulloblastoma  cell  lines  clustered  with  one  another,  but  failed  to  show 
similarity  to  primary  Group  3  samples  or  samples  from  any  other  subgroup. 

Super-enhancer  regulated  transcription  factors  implicate  Group  4  cellular  origins 

Among  subgroup-specific  SE  target  genes,  we  observed  an  enrichment  of  TFs  involved 
in  neuronal  development  (P~0.0001,  Fisher’s  exact  test;  Extended  Data  Figure  6a). 
-30%  of  these  SEs  are  previously  unreported  suggesting  a  lineage-specific  regulatory 
role.  Overall,  subgroup-specific  TFs  displayed  similar  patterns  of  expression,  enhancer 
motif  enrichment,  and  overlap  of  target  genes  (Extended  Data  Figures  6  &  7).  TFs  were 
also  enriched  in  subgroup-specific  SE  targets  as  compared  to  subgroup-specific  non-SE 
targets  (P~0.002,  Fisher’s  exact  test),  consistent  with  prior  observations  in  other  cancers 
that  SEs  regulate  key  TFs  required  for  tumour  identity  and  maintenance^^’^'*'^'*.  Given 
evidence  in  embryonic  stem  cells  that  pluripotency  master  regulator  TFs  (OCT4,  SOX2, 
NANOG)  are  driven  by  SEs  and  themselves  bind  to  and  establish  SEs^®,  we 
hypothesized  that  a  reverse  analysis  of  SEs  in  medulloblastoma  might  enable  a  de  novo 
reconstruction  of  tumour  identity-defining  TFs  and  their  associated  regulatory  circuitry, 
thereby  providing  novel  insights  into  medulloblastoma  origins. 

Pursuant  to  this  idea,  we  proposed  a  set  of  criteria  for  TF  inclusion  into  the  core 
regulatory  circuitry  of  medulloblastoma.  Specifically,  (i)  core  regulatory  circuitry  TFs  are 
SE-regulated  and  (ii)  the  TFs  themselves  bind  to  SEs  of  one  another  (Figure  5a).  For 
each  SE-regulated  TF,  these  criteria  can  be  quantified  through  a  measurement  of  the  in 
and  ouf  degree  of  regulation,  whereby  the  in  degree  represents  the  total  number  of  SE- 
regulated  TFs  that  bind  to  a  TF’s  SE,  and  the  out  degree  represents  the  total  number  of 
other  TF  SEs  bound  by  a  given  TF  (Figure  5a).  Using  these  criteria  in  the  poorly 
characterized  Group  4  medulloblastoma,  we  observe  interconnected  binding  at  the  SEs 
of  three  neuronal  TFs:  LMX1A,  LHX2,  and  EOMES  (Figure  5b).  Inspection  of  their 
respective  gene  loci  revealed  large  SEs  containing  clustered  binding  sites  for  these 
factors  present  only  in  Group  3  and  Group  4  (Figure  5b, c).  Additionally,  Group  4-specific 
enhancers  for  LMX1A,  LHX2,  and  EOMES  binding  sites  linked  those  TFs  with  Group  4- 
specific  target  genes  (Extended  Data  Figure  8).  Extending  regulatory  circuitry 
reconstruction  across  all  SE-associated  TFs  in  medulloblastoma,  we  identified  regulatory 
cliques  of  TFs  with  similar  patterns  of  in/out  degree,  strong  interconnectivity  via  motif 
binding,  and  higher  likelihoods  of  pairwise  protein/protein  interaction  and  motif  co¬ 
occurrence  at  enhancers  (Figure  5d,  Extended  Data  Figure  9).  This  reconstruction 
creates  for  the  first  time  a  candidate  core  regulatory  circuitry  in  each  subgroup,  and 
implicates  specific  sets  of  TFs  in  establishing  medulloblastoma  subgroup  identity  (Figure 
5d). 


Cellular  origins  for  WNT  and  SHH  medulloblastomas  have  been  experimentally 
established  using  mouse  models  genetically  engineered  to  aberrantly  activate  the  WNT 
and  SHH  signaling  pathways,  respectively,  in  distinct  hindbrain  stem/progenitor  cells 
during  development^'^°.  The  origins  of  Group  3  and  Group  4  medulloblastoma,  however, 
are  unknown  and  yet  essential  to  define  as  these  tumours  account  for  -60%  of  all 
diagnoses,  lack  targeted  therapies,  and  are  frequently  associated  with  a  poor  clinical 
outcome  secondary  to  current  standard  of  care^. 

Cell  identity  is  most  essentially  defined  by  the  activity  of  master  regulator  TFs.  In 
reprograming  and  trans-differentiation  studies,  the  activity  of  these  TFs  (e.g.  OCT4, 
SOX2,  and  NANOG  in  embryonic  stem  cells,  Pu.1  in  B  cells,  ATOH1  in  sensory  hair 
cells,  and  BARHL1/BARHL2  in  spinal  commissural  neurons)  is  sufficient  to  induce 
transitions  between  cell  states^®'"^^.  As  such,  we  hypothesized  that  the  regulatory  SE 
regions  governing  endogenous  expression  of  candidate  master  TFs  and  embedded  in 
the  core  regulatory  circuitry  of  medulloblastoma  subgroups  might  inform  the  cellular 
origins  of  the  disease  via  their  cell  type-specific  activity.  Examination  of  the  expression  of 
Lmx1a,  Eomes,  and  Lhx2  in  the  developing  mouse  hindbrain  (e13.5)  using  the  Mouse 
Allen  Brain  Atlas  Database  showed  spatiotemporal  patterns  of  restricted  expression  in 
the  nuclear  transitory  zone  (NTZ),  an  assembly  point  for  immature  deep  cerebellar  nuclei 
(DCN)  (Figure  5e).  DCN  residing  in  the  NTZ  at  this  time  point  are  predominantly 
glutamatergic  projection  neurons  that  originate  from  earlier  progenitors  of  the  rhombic 
lip,  a  transient  structure  producing  progenitors  with  distinct  cellular  fates,  including  DCN 
and  cerebellar  granule  neurons"^^  (Figure  5g).  Immunohistochemical  staining  for  Lmx1a 
and  Eomes  at  the  same  developmental  stage  (e13.5)  recapitulated  these  findings 
(Figure  5f).  To  investigate  the  spatiotemporal  distribution  of  cells  in  which  Group  4- 
specific  TFs  are  regulated  by  their  SEs,  we  cloned  constituent  regions  of  the  Lmx1a  SE 
into  a  LacZ  reporter  construct,  introduced  the  reporter  into  the  e11.5  developing  mouse 
hindbrain  by  exo  utero  electroporation,  and  assayed  enhancer  activity  via  x-gal  staining 
at  48  hours  post-transfection  (Figure  5g).  X-gal  staining  revealed  spatially  restricted 
activity  of  the  Lmx1a  SE  reporter  in  the  developing  cerebellum  (Figure  5h).  These 
findings  validate  the  specific  activity  of  the  Lmx1a  SE  observed  in  Group  4 
medulloblastoma  and  implicate  precursors  of  glutamatergic  DCN  as  potential  cells-of- 
origin  for  this  subgroup.  Finally,  these  findings  establish  SE  core  regulatory  circuitry  as  a 
novel  method  to  infer  cell  of  origin  for  poorly  classified  primary  tumours. 

Discussion  (450  words): 

We  describe  the  active  medulloblastoma  enhancer  landscape  across  a  series  of  28 
fresh-frozen,  treatment-naive  tissue  samples  and  three  cultured  cell  lines,  to  our 
knowledge  representing  the  largest  such  dataset  for  any  single  cancer  entity.  Our  data 
reveal  dramatic  divergence  between  primary  tumour  and  tumour  cell  line  material  and 
uncover  considerable  c/s-regulatory  element  heterogeneity  between  subgroups  of  the 
disease  that  would  be  overlooked  and  unsubstantiated  in  series  limited  to  just  a  few 
cases. 

Clinically  relevant  medulloblastoma  subgroups  are  principally  defined  based  on  their 
underlying  transcriptional  profiles.  Differentially  regulated  medulloblastoma  enhancers 
and  SEs  are  here  shown  to  recapitulate  these  subgroups,  and  importantly  extend  our 
understanding  of  this  disease  to  inferences  regarding  cell  specification  and  actionable 
tumour  dependencies.  Biological  themes  and  signaling  networks  extracted  from 
transcriptional  data  have  served  as  the  primary  source  of  annotation  for 


medulloblastoma  subgroups,  with  WNT  and  SHH  subgroups  characterized  by  activation 
of  their  respective  signaling  pathways,  and  Group  3  and  Group  4  recognized  for  their 
GABAergic  and  glutamatergic  expression  phenotypes,  respectively.  Although  these  data 
provide  a  functional  and  phenotypic  annotation  of  medulloblastoma,  they  fail  to  articulate 
the  cell  of  origin  and  developmental  identity  of  individual  subgroups.  Using  a  reverse 
analysis  of  the  medulloblastoma  chromatin  landscape  starting  at  the  level  of 
differentially-regulated  enhancers  and  SEs,  we  have  reconstructed  a  model  of  the  core 
regulatory  circuitry  inherent  to  medulloblastoma  subgroups,  and  inferred  master 
transcriptional  regulators  responsible  for  subgroup-specific  divergence.  The  majority  of 
these  master  regulator  TFs  were  not  previously  implicated  in  medulloblastoma 
developmental  biology,  nor  were  they  visible  amongst  transcriptionally-derived  gene  sets 
dominated  by  aberrant  signaling  and  overwhelming  phenotypic  signatures.  Through 
tracing  the  spatiotemporal  activity  of  a  subset  of  Group  4  master  TFs,  these  studies 
identified  DCN  of  the  cerebellar  NTZ,  or  plausibly  their  earlier  precursors  originating  from 
the  rhombic  lip,  as  putative  cells-of-origin  for  this  large  subgroup  of  patients.  Together 
these  approaches  establish  a  framework  for  the  inference  of  tumour  cell  of  origin  through 
enhancer  core  regulatory  circuitry  mapping. 

Understanding  the  cellular  origins  of  cancer  has  broad  implications  for  the  understanding 
and  treatment  of  the  disease'*'*.  Numerous  cancers,  especially  those  of  the  immune 
compartment  are  treated  through  targeting  of  the  lineage  {e.g.  anti-B  cell  therapies)'*®  '*®. 
As  medulloblastoma  is  believed  to  originate  from  cell  populations  that  normally  exist 
ephemerally  during  development,  targeting  the  aberrant  persistence  of  tumour  cells  from 
these  lineages  may  represent  a  novel  therapeutic  strategy.  Moreover,  elucidation  of  core 
regulatory  circuitry  implicates  upstream  signaling  dependent  regulators  of  master  TFs, 
their  co-activators,  and  their  downstream  effectors  as  potential,  rational  subgroup- 
specific  therapeutic  targets.  These  insights  demonstrate  the  critical  importance  of 
epigenetic  analyses  of  primary  tumours  as  opposed  to  cell  line  model  systems  and 
highlight  the  broad  utility  of  core  regulatory  circuitry  inference  especially  in  poorly 
characterized  and  clinically  diverse  malignancies. 

Methods  Summary 

All  patient  material  included  in  this  study  was  collected  after  receiving  informed  consent 
from  the  patients  and  their  families.  Medulloblastoma  samples  were  collected  at  first 
resection,  before  adjuvant  chemotherapy  or  radiotherapy.  Subgroup  assignments  were 
made  using  the  lllumina  450K  DNA  methylation  array  as  described'*^.  Chromatin 
extraction  and  library  preparation  for  ChIP-seq  of  H3K27ac  and  BRD4  were  performed 
at  ActiveMotif  (Carlsbad,  CA)  using  proprietary  protocols.  Alignment  and  filtering  of 
ChIP-seq  data  was  performed  as  described®.  H3K27ac  enhancer  peaks  were  called 
using  MACS^®.  H3K27ac  peaks  were  classified  as  being  subgroup-specific  or  as 
common  enhancers  by  first  calculating  H3K27ac  enrichment  on  the  merged  peaks 
followed  by  ANOVA  and  k-means  clustering.  Target  gene  identification  of  enhancers 
was  performed  by  correlating  H3K27ac  enrichment  at  the  enhancers  with  expression 
levels  of  genes  located  in  the  same  TAD®*  as  the  enhancers.  Candidate  gene(s) 
showing  the  highest  correlation  were  selected  as  the  putative  target(s)  of  the  enhancer. 
Gene  Ontology/Pathway  analysis  of  enhancer-gene  targets  was  performed  using  the 
ClueGO  plugin  for  cytoscape'*®.  SEs  were  called  using  ROSE2 
(https://qithub.com/BradnerLab/pipeline)'*^  and  subgroup  specificity  of  super-enhancers 
were  assigned  via  ranking  average  H3K27ac  signal  across  the  subgroups. 
Medulloblastoma  core  regulatory  circuitry  analysis  was  performed  by  calculating  inward 


and  outward  degree  regulation  of  SE-regulated  TFs.  Reporter  assays  for  validating 
enhancers  in  vivo  were  performed  by  exo  utero  electroporation  of  reporter  constructs 
into  the  hindbrain  of  murine  CD1  embryos  (e11.5).  Enhancer  activity  was  measured  by 
x-gal  staining  of  transfected  cerebella.  Endogenous  expression  of  candidate  TFs  was 
determined  by  querying  the  Allen  Brain  Atlas  Data  Portal  (http://developinqmouse.brain- 
map.orq)  or  by  immunohistochemistry  performed  on  the  murine  embryonic  hindbrain 
(e13.5). 
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Figure  Legends: 

Figure  1.  The  enhancer  landscape  of  primary  medulloblastoma. 

(a)  Experimental  workflow  for  studying  enhancers  and  super-enhancers  in  primary 
medulloblastomas. 

(b)  H3K27ac  ChIP-seq  data  across  all  28  primary  medulloblastoma  samples  from 
our  series  showing  a  highly  active  enhancer  at  the  OTX2  locus,  especially  in 
Group  3  and  Group  4. 

(c)  Scatter  plot  comparing  the  enrichment  (log2)  of  H3K27ac  signal  versus  BRD4 
signal  at  medulloblastoma  enhancers  (n=78,516)  as  determined  by  ChIP-seq. 

(d)  Scatter  plot  comparing  the  enrichment  (log2)  of  H3K27ac  signal  versus  DNA 
methylation  at  medulloblastoma  enhancers  (n=78,516)  as  determined  by  ChIP- 
seq  and  WGBS,  respectively. 

(e)  RNA-seq  data  showing  Group  3-specific  enhancer  RNA  (eRNA)  expression 
(lower  left)  overlapping  a  Group  3-specific  /WYC  enhancer  (upper  left)  in  a  subset 
of  Group  3  and  Group  4  medulloblastomas.  MYC  gene  expression  (RPKM)  is 
also  shown  for  the  same  subset  of  cases  (lower  right). 

(f)  Venn  diagram  showing  the  overlap  of  medulloblastoma  enhancers  with  those 
reported  by  ENCCDE  and  Roadmap. 

(g)  Pie  chart  summarizing  the  results  presented  in  (f). 

Figure  2.  Differentially  regulated  enhancers  and  enhancer-gene  assignments  in 
medulloblastoma  subgroups. 

(a)  ANCVA  classification  of  medulloblastoma  enhancers  displayed  as  a  pie  chart. 

(b)  Pie  chart  showing  the  distribution  of  differentially-regulated  enhancers  among 
medulloblastoma  enhancer  classes. 

(c)  K-means  clustering  of  differentially-regulated  medulloblastoma  enhancers 
(n=20,406). 

(d)  Bar  plot  showing  the  proportion  of  enhancer/gene  assignments  to  N  enhancers. 

(e)  Bar  plot  displaying  the  proportion  of  enhancer/gene  assignments  to  N  genes. 


(f)  Bar  plot  summarizing  the  proportion  of  enhancer-gene  assignments  to  N 
subgroups. 

(g)  WNT  (E2)  and  SHH  (E3)  subgroup-specific  enhancers  inferred  to  regulate  ALK. 

(h,  i)  Scatter  plots  correlating  sample-matched  gene  expression  (RPKM,  x-axis)  of 

ALK  with  H3K27ac  enrichment  (log2:  y-axis)  for  the  WNT-specific  (E2)  and  SHH- 
specific  (E3)  enhancers  shown  in  (g). 

Figure  3.  Functional  characterization  of  enhancer/gene  targets  in 
medulloblastoma  subgroups. 

(a)  Waterfall  plot  discriminating  the  top  1,000  Group  3  and  Group  4  subgroup- 
specific  enhancers  as  defined  by  total  H3K27ac  signal.  The  distribution  of 
assigned  targets  in  Group  3,  Group  4,  and  shared  Group  3-4  targets  are  shown 
below  the  waterfall. 

(b,  c)  Functional  annotation  of  target  genes  assigned  to  Group  3  and  Group  4 
subgroup-specific  enhancers  based  on  their  significant  overlap  with  gene  sets 
annotated  in  Gene  Ontology  (GO  Biological  Process)  and  pathway  databases 
(KEGG,  Reactome). 

(d)  Convergence  of  Group  3-specific  enhancers  on  TGFp  pathway  genes. 
Subgroup-specific  enhancers  are  summarized  as  nodes  according  to  their 
respective  medulloblastoma  enhancer  class  -  Group  3,  Group  4,  and  shared 
Group  3-4  -  with  edges  representing  individual  enhancer/TGFp  pathway  gene 
assignments. 

Figure  4.  Medulloblastoma  super-enhancers  define  subgroup-specific  identity. 

(a)  Unsupervised  hierarchical  clustering  of  primary  medulloblastomas  and  cell  lines 
using  H3K27ac  signal  calculated  at  all  SEs  identified  in  each  individual  sample. 

(b)  Meta  tracks  of  H3K27ac  ChIP-seq  signal  for  the  conserved  ZIC1/ZIC4  SE  locus. 
Expression  (mean  RPKM)  for  both  ZIC4  (left)  and  ZIC1  (right)  is  displayed  as  bar 
graphs  to  the  right  of  each  H3K27ac  track. 

(c)  Line  plot  showing  the  enhancer  rank  for  the  conserved  ZIC1/ZIC4  SE  locus 
across  all  samples  according  to  subgroup. 

(d)  Ranked  plots  of  enhancers  defined  across  composite  H3K27ac  landscapes  of 
WNT,  SHH,  Group  3,  and  Group  4  medulloblastomas.  Enhancers  are  ranked  by 
increasing  group  average  H3K27ac  signal  (rpm).  The  cut-off  discriminating 
typical  enhancers  (TEs)  from  super-enhancers  (SEs)  is  shown  as  a  dashed  line. 
Select  genes  associated  with  SEs  in  each  subgroup  are  highlighted  and  shaded 
according  to  enhancer  class  specificity. 

(e)  Line  plots  showing  the  enhancer  rank  for  candidate  SE  loci  across  all  samples 
according  to  subgroup.  Examples  of  subgroup-specific  (WNT=A/K'D2, 
SHH=PCNT,  Group  3=HLX,  Group  4=SNCAIP)  SEs  are  shown. 

(f)  Meta  tracks  of  H3K27ac  ChIP-Seq  signal  across  medulloblastoma  subgroups  for 
the  loci  described  in  (e).  The  y-axis  shows  ChIP-Seq  signal  (rpm/bp)  for  each 
individual  sample  (shaded  regions)  with  the  average  signal  across  the  group 
shown  in  a  line.  The  x-axis  depicts  genomic  position  with  SE  boundaries 
demarcated  as  rectangles.  Bar  graphs  shown  to  the  right  of  each  H3K27ac  track 
summarize  the  expression  (mean  RPKM)  of  the  relevant  candidate  genes  as 
determined  by  RNA-seq. 


Figure  5.  Super-enhancers  define  medulloblastoma  regulatory  circuitry  and 
identify  putative  cellular  origins. 


(a)  Methodology  for  inferring  transcriptional  regulatory  circuitry  driven  by 
medulloblastoma  SEs. 

(b)  Subset  of  the  Group  4-specific  transcriptional  network  predicted  to  be  driven  by 
LMX1A,  LHX2,  and  EOMES. 

(c)  H3K27ac  ChIP-seq  meta  tracks  for  the  SE-regulated  TFs  LMX1A,  LHX2,  and 
EOMES.  Locations  of  enriched  motifs  for  each  of  the  respective  TFs  are 
highlighted  at  the  top  of  the  panel. 

(d)  Subgroup-specific  TF  circuitry.  Nodes  are  TFs  associated  with  a  SE  in  a 
subgroup-specific  context.  Edges  indicate  co-regulating  TFs  as  defined  by 
enrichment  of  TF  binding  motifs  in  respective  regulatory  regions. 

(e)  In  situ  hybridization  data  showing  highly  localized  expression  of  Lmx1a  (upper 
panel)  and  Eomes  (middle  panel)  in  the  embryonic  cerebellum  at  e13.5.  Red 
arrows  indicate  highly  specific  focal  expression  of  both  TFs  in  deep  cerebellar 
nuclei  (DCN)  of  the  nuclear  transitory  zone  (NTZ).  Atoh1  expression  (lower 
panel)  is  shown  at  the  same  developmental  time-point  to  serve  as  a  marker  of 
the  external  granule  layer  (EGL). 

(f)  Immunofluorescence  microscopy  for  Lmx1a  (upper  panel)  and  Eomes  (middle 
panel)  performed  on  sagittal  sections  of  the  e13.5  murine  cerebellum.  Red 
arrows  indicate  highly  localized  expression  of  both  TFs  in  DCN  of  the  NTZ. 
Atoh1  staining  (lower  panel)  is  shown  at  the  same  developmental  time-point  to 
serve  as  a  marker  of  the  EGL.  Neighbouring  sections. 

(g)  Strategy  for  validating  medulloblastoma  subgroup-specific  enhancers  in  vivo. 
An  Atlas  (sagittal)  of  the  e13.5  murine  cerebellum,  highlighting  the  main  cell 
types  and  compartments  contributing  to  the  cerebellar  aniage  at  this 
developmental  stage  is  included  as  a  reference. 

(h)  X-gal  staining  of  an  embryonic  hindbrain  (e13.5)  transfected  with  the  reporter 
construct  shown  in  (g)  containing  constituents  of  the  Group  4-specific  Lmx1a 
SE.  Red  arrow  indicates  cells  positive  for  in  vivo  enhancer  activity  driving  LacZ 
expression.  Dorsal  view. 
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Figure  1.  The  enhancer  landscape  of  primary  medulloblastoma 
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Figure  2.  Differentially  regulated  enhancers  and  enhancer-gene  assignments  in 
medulloblastoma  subgroups. 
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Figure  3.  Functional  characterization  of  enhancer-gene  targets  in 
medulloblastoma  subgroups. 
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Figure  4.  Medulloblastoma  super-enhancers  define  subgroup-specific  identity 
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Figure  5.  Super-enhancers  define  meduilobiastoma  reguiatory  circuitry  and 
identify  putative  ceilular  origins. 
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Myc  deregulation  is  a  hallmark  oncogenic  event  where 
overexpression  of  the  transcription  factor  gives  rise  to 
numerous  tumorigenic  phenotypes.  The  complex  con¬ 
sequences  of  Myc  deregulation  have  prevented  clear 
mechanistic  interpretations  of  its  function.  A  synthesis 
of  recent  experimental  observations  offers  a  consensus 
on  the  direct  transcriptional  function  of  Myc:  when 
overexpressed,  Myc  broadly  engages  the  established 
euchromatic  c/s-regulatory  landscape  of  the  cell,  where 
the  factor  generally  amplifies  transcription.  The  level  of 
Myc  binding  at  target  genes  and  the  transcriptional 
output  are  differentially  modulated  by  additional  regu¬ 
lators,  including  Miz1.  Targeting  Myc  oncogenic  activity 
will  require  an  understanding  of  whether  amplification 
promotes  tumorigenesis  and  the  consequences  of  am¬ 
plification  in  tumors  adapted  to  oncogenic  Myc. 

Myc  transcription  in  cancer 

The  intense  interest  in  Myc  proteins  stems  from  their 
pervasive  role  in  the  genesis  of  human  tumors.  A  large 
hody  of  evidence  has  established  that  level  of  one  out  of 
three  Myc  proteins  (c-Myc,  N-Myc,  or  L-Myc)  is  enhanced 
and  its  expression  released  from  its  normally  tight  depen¬ 
dence  on  growth  factors  in  a  large  fraction  of  all  human 
tumors  [1,2].  The  three  Myc  genes  are  differentially 
expressed  during  development,  but  the  proteins  are  func¬ 
tionally  equivalent  in  most  biological  systems,  allowing  us 
to  focus  on  c-Myc  (called  ‘Myc’  from  here  on)  in  this  review 

[3]  .  Multiple  experiments  in  transgenic  models  of  human 
tumors  show  that  this  ‘deregulated’  expression  of  Myc 
proteins  promotes  tumorigenesis  and  that  tumors  gener¬ 
ated  by  different  oncogenes  depend  on  elevated  Myc  levels 

[4] .  Accordingly,  several  proof-of-principle  studies  suggest 
that  targeting  Myc  proteins  have  considerable  benefit  for 
tumor  therapy  [5,6].  This  hope  awaits  clinical  confirma¬ 
tion,  because  therapies  directed  against  Myc  proteins 
have  yet  to  enter  the  clinic.  Myc  proteins  are  nuclear 
proteins  and,  as  such,  may  have  several  functions;  howev¬ 
er,  there  is  consensus  in  the  field  that  one  major  activity  of 
Myc  proteins  is  to  regulate  transcription  [7].  This  has 
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stimulated  a  vigorous  debate  about  why  Myc-dependent 
transcription  is  so  important  for  tumorigenesis.  Here,  we 
review  the  current  state  of  this  debate. 

The  evolution  of  Myc 

One  of  the  best-understood  Myc  transcription  factor  net¬ 
works  is  that  of  Drosophila.  dMyc,  a  basic  helix-loop-helix 
leucine  zipper  protein  (bHLH-Zip),  binds  together  with  a 
another  bHLH-Zip  protein,  dMax,  to  target  genes  via  a 
conserved  sequence  element  termed  an  ‘E-box’  (CACGTG) 
and  to  activate  transcription  [8].  The  dMyc  antagonist 
dMnt,  also  pairs  with  dMax  to  repress  transcription  from 
the  same  DNA  sequence  [9].  Deletion  of  dMnt  largely 
rescues  the  developmental  defects  of  dMyc  deletion,  pro¬ 
viding  strong  evidence  that  transcriptional  activation  is 
the  critical  biological  function  of  dMyc  [10].  Most  target 
genes  of  dMyc  and  dMax  encode  proteins  involved  in  ENA 
and  protein  synthesis  [11].  Consistently,  dMyc  proteins 
stimulate  cell  growth,  but  not  cell  proliferation,  and  loss-of- 
function  alleles  show  a  ‘minute’  phenotype  [8,12].  Transfer¬ 
ring  this  model  to  mammalian  tumorigenesis  would 
suggest  that  the  oncogenic  activity  of  Myc  is  due  to  the 
deregulated  and  constitutive  activation  of  its  target  genes. 
According  to  this  model,  sustained  high  levels  of  Myc  in 
human  tumors  over-drive  RNA  biogenesis  and  protein 
translation,  enabling  growth  factor-independent  cell 
growth.  Multiple  observations  are  consistent  with  this 
model:  for  example,  Myc-driven  lymphomas,  in  contrast 
to  their  normal  counterparts,  rely  on  supraphysiological 
rates  of  protein  translation  to  support  growth  [13]. 

Similar  to  Drosophila  dMyc,  mammalian  Myc  proteins, 
when  dimerized  with  Max,  are  sequence-specific  DNA-bind- 
ing  proteins  that  bind  to  E-boxes  [7].  Unlike  prokaryotic 
repressors  that  discriminate  operator  from  nonoperator 
DNAs  by  approximately  10®“®-fold,  Max  homodimers  and, 
most  likely,  Myc/Max  heterodimers  have  considerable  affin¬ 
ity  for  variant  E-boxes  and  even  generic  DNA  [14].  There¬ 
fore,  it  is  almost  certain  that  Myc/Max  heterodimers  have 
non-negligible  affinity  for  most  promoters.  Reflecting  its 
continuum  of  association  to  targets,  the  intranuclear 
d3mamics  of  Myc  suggest  engagement  with  intranuclear 
partners,  indicating  that  Myc  is  globally  distributed,  is 
almost  always  bound  directly  or  indirectly  with  chromatin, 
and  is  available  to  all  active  genes  [15,16].  While  broadly 
anticipated,  there  is  currently  only  limited  formal  evidence 
that  transcriptional  activation  is  the  critical  oncogenic 


Trends  in  Cell  Biology  xx  (2014)  1-8  1 


ARTICLE  IN  PRESS 


s  1 100;  No.  of  Pages  8 


Review 


Trends  in  Cell  Biology  xxx  xxxx,  Vol.  xxx,  No.  x 


function  of  Myc.  One  example  is  the  finding  that  inactivation 
of  Mga,  a  repressive  partner  protein  of  Max,  occurs  com¬ 
monly,  but  mutually  exclusively  with  focal  amplification  of 
the  Myc  gene  in  lung  cancer,  suggesting  that  Mga  inactiva¬ 
tion  is  pathologically  equivalent  to  Myc  overexpression  in 
this  entity  [17]. 

Ques  to  how  Myc  regulates  gene  expression  are  provided 
in  the  list  of  co-activators  it  recruits  to  induce  transcription: 
a  conserved  sequence  in  Myc  (MycBoxII)  binds  the  adaptor 
protein  Trans&rmation/transcription  domain-associated 
protein  (TRRAP)  via  TRRAP  GCN  5-containing  histone  acet- 
ylase  complexes  [18,19].  Myc  also  recruits  the  p300  histone 
acetylase  and  other  histone-modifying  enaymes  to  its  target 
sites,  although  the  requisite  interaction  simfaces  on  Myc  are 
less  well  defined  [20,21],  Myc  also  interacts  with  TATA  box 
binding  protein  (TBP)  [22]  and  the  chromatin  regulators 
Tip60  [23]  and  switch/sucrose  nonfermentable  (SWI/SNF) 
[24],  and  inhibits  the  H3  trimethyl  K4  (H3K4Me3)  demethy- 
lase  dKDM5/LID  [25],  The  ability  of  Myc  to  interact  with 
and  recruit  a  multitude  of  general  chromatin  regulators 
is  consistent  with  its  global  role  in  r^ulating  chromatin 
structure  because  Myc  activity  has  been  broadly  linked  with 
large  nuclei/nucleoli  containing  high  levels  of  euchromatic 
and  hyperacetylated  chromatin  [18,26]. 

Myc  may  also  interact  directly  with  the  transcription 
apparatus.  Most  notably,  the  extreme  amino-terminus  of 
Myc  binds  (yclin  T1  in  vitro  [27].  Cyclin  T1  is  part  of  the 
positive  transcription  elongation  factor  b  (P-TEFb)  elonga¬ 
tion  complex,  which  phosphorylates  the  C-terminal  do¬ 
main  of  RNA  polymerase  at  serine  2  and  thereby 
promotes  the  transition  fi-om  promoter  binding  to  produc¬ 
tive  elongation.  The  relevance  of  this  finding  stems  fi"om 
observations  that  Myc  can  stimulate  transcription  at  a  step 
after  initiation  of  transcription  and  that  RNA  pol5anerase 
II  acciunulates  in  the  body  of  transcribed  genes  after 


activation  of  Myc  [28-30].  Fmthermore,  ectopic  expression 
of  Myc  globally  enhances  phosphorylation  of  RNA  poly¬ 
merase  [30,31].  Together,  these  observations  suggest  that 
the  release  of  RNA  polymerase  fi-om  a  paused  position  is  a 
key  mechanism  of  transcriptional  activation  by  Myc  [30] 
(Figure  1).  There  is  no  doubt  that  other  effectors  of  Myc 
remain  to  be  defined.  How  Myc  exploits  these  partners 
individually  or  in  combination  at  any  given  promoter  is  not 
yet  known. 

The  general  amplifier  model 

The  identification  of  Myc  as  a  transcriptional  activator 
raised  the  expectation  that  enumeration  of  the  direct 
target  genes  of  Myc  would  provide  a  list  of  critical  down¬ 
stream  targets  and  biological  processes  that  mediate  the 
physiological  and  oncogenic  functions  of  Myc.  This  expec¬ 
tation  prompted  a  series  of  studies  to  identify  the  Myc- 
regulated  genes  by  comparing  RNA  expression  profiles  and 
the  genome-wide  map  of  Myc-boimd  chromatin  using 
microarray  or  next  generation  sequencing  technology. 
Interrogation  of  a  range  of  cell  t3T)es  surprisingly  revealed 
mammalian  Myc  fjroteins  at  nearly  aU  promoters  in  open 
chromatin  [32].  This  was  not  only  true  for  genes  tran¬ 
scribed  by  RNA  polymerase  II,  but  also  for  genes  tran¬ 
scribed  by  RNA  polymerases  I  and  HI  (although  rRNA 
promoters  tend  to  bind  less  Myc  compared  with  other 
promoters)  [33,34].  For  RNA  polymerase  Il-transcribed 
genes,  the  presence  of  histone  marks  indicative  of  open 
chromatin,  in  particular  H3K4Me3,  highly  correlated  with 
Myc  binding  [35].  When  overexpressed,  Myc  proteins  also 
bound  to  many  enhancers,  again  correlating  with  histone 
marks  indicative  of  activity.  Given  that  many  genes  encod¬ 
ing  proteins  involved  in  translation  (e.g.,  ribosomal  protein 
genes)  are  highly  transcribed  fium  promoters  in  CpG 
islands  where  E-boxes  are  more  fii^quently  found  in  an 
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Figure  1 .  Schematic  of  transcription  elongation  promoting  activity  of  MYC.  M YC/MAX  bind  to  E-box  (CACGTG)  sequences  at  enhancers  and  promoters  and  recruit  the 
elongation  factor  complex  positive  transcription  elongation  factor  b  (P-TEFb),  whk:h  promotes  elongation  in  part  by  phosphorylating  the  C-terminal  domain  of  RNA 
Polymerase  IMRNA  Pol  II). 
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open  chromatin  conformation,  these  results  parallel  those 
obtained  in  Drosophila. 

Collectively,  these  data  have  led  to  an  exciting  model  in 
^vhich  Myc  proteins  act  as  general  amplifiers  that  globally 
enhance  transcription;  the  specificity  and  degree  of  ampli¬ 
fication  at  individual  genes  is  determined  by  the  exact 
pattern  of  open  promoters  that  pre-exists  in  a  cell  even 
before  Myc  activation  [36,37].  In  this  model,  the  unique 
biological  role  of  Myc  stems  from  the  ability  to  accelerate 
transcriptional  activation  (i.e.,  releasing  a  paused  RNA 
polymerase  from  the  promoter)  rather  than  by  specif3dng  a 
particular  program  of  gene  expression  [30].  Strong  support 
for  this  model  comes  from  the  analysis  of  primary  B  cells 
[36].  Upon  stimulation  of  resting  B  cells,  Myc  binds  to  all 
promoters  of  transcriptionally  active  genes  in  these  cells 
and  Myc  binding  follows  the  pattern  of  loading  of  RNA 
polymerase  II  that  pre-exists  in  resting  low  Myc  cells.  In 
response  to  mitogen  stimulation,  there  is  a  global  increase 
in  mRNA  levels  that  depends  on  Myc,  which  has  been 
confirmed  by  examining  the  mRNA  levels  of  multiple 
individual  genes.  There  are  indications  that  the  ‘global 
amplifier’  model  operates  not  only  in  primary  cells,  but  also 
in  tumors.  Comparing  normal  B  cells  with  Myc-driven  B 
cell  lymphomas  revealed  a  strong  increase  in  total  cellular 
mRNA  content  in  lymphomas  [38] .  ChIP  sequencing  shows 
that  Myc  binds  to  90%  of  all  open  promoters  in  lymphomas; 
an  easy  explanation  might  be  that  Myc  broadly  enhances 
expression  of  the  genes  to  which  it  binds  and  this  increases 
cellular  mRNA  content.  Similarly,  the  comparison  of  hu¬ 
man  tumor  cell  lines  with  different  Myc  levels  suggests 
that  Myc-driven  amplification  accounts  for  differences  in 
mRNA  content  between  different  tumor  cells  [37]. 

According  to  a  strict  amplifier  model,  the  direct  effect  of 
Myc  on  promoters  is  always  positive.  The  effect  of  Myc  is 
not  equally  strong  for  each  gene  because  promoters  and 
relevant  enhancers  differ  in  their  affinity  for  Myc:  promot¬ 
er  binding  dictates  the  differences  in  response.  Moreover, 
the  overall  promoter  response  to  Myc  is  not  directly  pro¬ 
portional  to  the  amount  Myc  recruited  there.  Furthermore, 
at  high  and  oncogenic  levels,  Myc  ‘invades’  enhancers  and 
promoters,  enabling  the  pathological  regulation  of  genes 
that  bind  little  Myc  at  physiological  levels  [37-39].  There¬ 
fore,  the  general  amplifier  model  can  account  for  the 
observation  that  oncogenic  levels  of  Myc  regulate  a  specific 
set  of  genes  in  particular  tumors  [40] .  The  model  can  also 
account  for  apparent  Myc-dependent  repression.  Specifi¬ 
cally,  if  expression  of  most  mRNAs  is  amplified  and,  there¬ 
fore,  if  total  mRNA  content  increases  in  a  cell  upon  Myc 
stimulation,  then  the  relative  amount  of  mRNAs  tran¬ 
scribed  from  genes  that  have  low-affinity  promoters  and, 
thus,  are  poorly  activated  genes  will  decrease  and,  depend¬ 
ing  on  the  mode  of  normalization,  these  genes  will  appear 
repressed  [41].  Hence,  whether  the  general  amplifier  mod¬ 
el  can  account  for  Myc-dependent  repression  is  intimately 
tied  to  the  question  of  whether  Myc  induces  a  general 
increase  in  cellular  RNA  levels  and,  by  inference,  whether 
cells  grow  as  a  resnlt  of  Myc  activation:  if  cells  grow, 
seemingly  repressed  genes  may  in  fact  be  low-affinity 
Myc  targets  that  are  relatively  ‘left  behind’  [41] .  Myc  drives 
both  growth  and  proliferation  of  mammalian  cells,  unlike 
in  Drosophila  [42].  Given  that  cell  proliferation  and 


division  reduce  cell  size,  some  of  the  apparent  Myc-depen¬ 
dent  repression  may  reflect  the  postmitotic  partition  of 
cellular  contents  compared  with  nondividing  cells;  there¬ 
fore,  a  complete  description  of  amplification  will  be  sensi¬ 
tive  to  cell  cycle  parameters.  Furthermore,  if  Myc  acts  at 
transcription,  then  gauging  Myc  action  from  mRNA  abun¬ 
dance  is  imprecise  because  mRNA  synthesis  becomes  con¬ 
flated  with  mRNA  half-lives,  and  cell  cycle  parameters. 
Yet,  this  too  becomes  even  more  complicated  because  the 
factors  that  regulate  mRNA  stability  and  the  cell  cycle  may 
themselves  be  important  direct  or  indirect  Myc  targets. 

An  important  aspect  of  the  amplifier  model  is  the  pos¬ 
tulate  that,  in  any  given  cell  type,  all  transcription  initiat¬ 
ing  genes  are  potential  targets  of  Myc:  individual 
promoters  differ  quantitatively  in  their  affinity  for  Myc, 
but  no  sharp  cut-off  discriminates  Myc  targets  from  non¬ 
targets.  Occupancy  at  each  promoter  changes  as  Myc  con¬ 
centrations  vary  and  this  change  defines  the  degree  of  the 
transcriptional  response.  We  suggest  that  there  is  a  con¬ 
tinuum  of  activity  as  Myc  flickers  on  and  off  of  weakly 
bound,  weakly  expressed  promoters,  but  stays  longer  or 
more  frequently  at  high  output  promoters  (Figure  2).  This 
view  of  Myc  as  a  general  regulator  capable  of  interacting 
with  all  active  promoters  has  implications  for  the  way  that 
Chip  seqnencing  data  need  to  be  interpreted,  because  it 
suggests  that  cut-offs  that  are  introduced  by  analysis  (e.g., 
peak-calling)  programs  generate  artificial  distinctions  be¬ 
tween  target  and  nontarget  genes. 

This  quantitative  view  also  suggests  that,  in  some 
experimental  settings,  the  response  of  some  genes  to 
Myc  plateaus  as  occupancy  by  Myc  at  promoters  and 
enhancers  and  corresponding  transcriptional  output  is 
saturated.  Whereas  Myc  abundance  is  low  in  most,  but 
not  all  primary  cells  (luminal  cells  in  mammary  ducts  are  a 
notable  exception),  recent  estimates  for  Myc  in  tumor  cells 
appear  high  enough  to  allow  the  saturation  of  some  pro¬ 
moters  [37].  A  recent  study  provides  evidence  that  high- 
affinity  promoters  are  Myc  saturated  in  proliferating  cells 
and,  therefore,  that  further  increases  in  Myc  only  increase 
occupancy  and  gene  expression  at  low-affinity  promoters 
[39].  As  a  result,  gene  expression  patterns  that  discrimi¬ 
nate  cells  with  oncogenic  Myc  levels  from  normally  prolif¬ 
erating  cells  are  different  from  the  gene  expression 
changes  that  occur  as  Myc-levels  increase  during  the 
physiological  transition  from  quiescence  to  proliferation. 

These  considerations  allow  the  conclusion  that,  even 
within  the  framework  of  the  general  amplifier  model,  the 
observed  effects  of  Myc  on  gene  expression  depend  on:  (i) 
the  precise  Myc  levels  that  are  being  compared;  (ii)  wheth¬ 
er  cells  grow  upon  Myc  activation;  and  (iii)  the  precise 
pattern  of  ongoing  transcription  that  determines  which 
genes  are  accessible  to  Myc. 

Direct  and  indirect  transcriptional  repression 

In  many  studies,  almost  as  many  genes  are  downregulated 
by  Myc  as  are  upregulated,  evoking  the  notion  that  Myc 
can  act  as  a  repressor  [38,43] .  Setting  aside  pseudodown- 
regulation  arising  from  normalization  issnes  (above),  such 
downregulation  can  stem  from  a  direct  and  context-depen¬ 
dent  repressive  function  of  Myc  that  is  unmasked  by 
factors  bound  locally  at  the  regulatory  sequences  of 
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Rgur«2.  Schematic  of  MYC  genomic  binding  at  (Ai  low,  (B|  n>edium,  and  (C|  high  concentrations  of  the  factor.  C/a-regulatory  elements  (boxes)  are  shaded  by  their  affinity 
for  Myc.  Myc  transcriptional  activity  at  genes  is  depicted  as  a  circle  sized  by  the  impact  of  Myc. 


particular  target  genes.  For  example,  Myc  blunts  the 
induction  of  the  cychn-dependent  kinase  4  inhibitor  B 
(CDKN2B)  and  CDKNIA  genes  in  response  to  transform¬ 
ing  growth  factor  (TGF)-3  signaling  by  a  direct  interaction 
with  the  Mizl  and  SmadS  proteins  at  both  promoters 
[44,45].  Alternatively,  repression  can  be  indirect;  as  a 
universal  amplifier,  Myc  will  increase  the  expression  of 
repressive  transcriptional  and  chromatin  components.  As 
repressors  are  upregulated,  downregulation  of  sets  of  tar¬ 
gets  may  ensue  as  the  system  seeks  a  new  steady  state.  For 
example,  Myc-upregulated  phosphatase  and  tensin  homo¬ 
log  (PTEN)  inactivates  Akt,  which  disinhibits  enhancer  of 
zeste  homolc^  2  (Ezh2),  promoting  widespread  polycomb 
repressive  complex  2  (PRC2)-mediated  repression  [46]. 
Furthermore,  miRNAs,  such  as  miR17-92  that  represses 
SIN3  transcription  regulator  family  member  B  (Sin3h), 
HMG-box  transcription  factor  1  (Hhpl),  suppressor  of 
vari^ation  4-20  homolc^  1  (Suv420hl),  and  B  cell  trans¬ 
location  gene  1  (Btgl)  to  regulate  chromatin  structure  as 
weU  as  the  apoptosis  facilitator  Bim  [47],  are  repressive 
mediators  of  Myc  action.  Direct  and  indirect  mechanisms 
may  become  entangled  as  feedforward  and  feedback  arms 
of  genetic  circuits  adjust  to  new  conditions.  In  principle, 
the  dynamics  of  the  temporal  profile  of  the  transcriptome 
could  help  to  discriminate  direct  fium  indirect  regulation, 
because  changes  requiring  newly  S5mthesized  protein 
intermediaries  are  expected  to  he  delayed  relative  to  the 
induction  of  Myc  protein.  However,  few  studies  have  had 
the  temporal  resolution  to  discriminate  direct  versus  indi¬ 
rect  targets.  Given  that  they  are  translation  independent, 
miRNAs  and  other  noncoding  RNAs  are  likely  to  act  after 
only  a  brief  delay.  This  rapid  onset  of  activity  is  likely  to 


hinder  the  discrimination  of  the  direct  and  indirect  effects 
of  Myc. 

Shaping  transcriptional  amplification 

Transcriptional  amplification  and  cell  growth  induced  by 
Myc  cannot  continue  indefinitely  and  unoppKjsed.  There¬ 
fore,  mechanisms  must  exist  that  adjust  transcriptional 
amplification  to  available  metabohc  resources  and  that 
provide  feedback  from  the  physiological  status  of  a  cell 
to  Myc  activity.  Consequently,  supraphysiological  levels  of 
Myc  would  be  expected  to  elicit  either  compensatory  mea- 
siues  or,  if  they  fail,  promote  cell  death.  As  Myc  saturates 
its  strongest  and  most  highly  expressed  sites  (these  being 
most  likely  to  be  physiological  targets),  it  may  patholog¬ 
ically  spin  over  onto  other  promoters  and  enhancers,  in¬ 
creasing  their  activity  and  overall  cellular  RNA 
production,  provided  that  the  cell  can  marshal  sufficient 
resoim:es  to  augment  net  transcription. 

There  is  evidence  for  feedback  mechanisms  that  limit 
Myc-dependent  transactivation  operating  at  both  physio¬ 
logical  and  supraphysiological  levels.  The  apoptosis 
elicited  by  overexpression  provides  an  example  of  an 
upper-boimd  limit  for  Myc  levels  [48].  As  another  example, 
the  ribosomal  protein  Lll,  encoded  by  a  direct  target 
gene  of  Myc,  inhibits  Myc-dependent  transactivation  and 
expression  [49,50].  Supraphysiological  levels  of  Myc  induce 
expression  of  the  Arf  tiunor  suppressor  protein,  which  in 
turn  interacts  with  Myc  and  inhibits  Myc-dependent  trans¬ 
activation  [51,52].  One  broadly  acting  mechanism  that 
limits  transcriptome  amplification  driven  by  supraphysio¬ 
logical  Myc  is  the  recruitment  of  the  zinc  finger  protein 
Mizl  at  many  promoters  [53].  In  primary  cells,  and  in  the 
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absence  of  high  levels  of  Myc,  Mizl  is  a  highly  sequence- 
specific  transcription  factor  that  activates  a  relatively 
small  number  of  targets  that  share  an  extended  consensus 
sequence  and  that  control  autophagy  [54] .  As  Myc  levels  rise 
to  supraphysiological  levels,  Mizl  joins  Myc  via  a  protein- 
protein  interaction  at  a  large  fi"action  of  promoters  that  then 
decreases  output  [39].  Mizl  is  imstable  and  continuously 
degraded  by  the  Huwel  ubiquitm  Ugase  [55];  upon  inhibi¬ 
tion  of  Huwel,  Mizl  broadly  acciunulates  at  Myc-bound 
promoters,  blunts  activation,  and  enhances  repression  by 
Myc  [55,56].  Interestingly,  the  Arf  protein  binds  to,  and 
inhibits,  Huwel  and  promotes  association  of  Myc  with 
Mizl  [57,58],  strongly  suggesting  that  Arf  and  Mizl  are 
part  of  a  common  stress  response  to  oncogenic  Myc  levels. 

A  second  mechanism  that  is  likely  to  limit  and  shape 
Myc  amplification  is  competition  between  Myc  and  other 
bHLH-Zip  dimers  for  binding  at  El-boxes  and  perhaps  at 
other  sites  (Figure  3).  Myc-Max  operates  against  Max 
homodimers  as  well  as  heterodimers  of  Max  with  repres¬ 
sive  members  of  the  bHLH-Zip  family  (Mxdl-4,  Mnt,  etc) 
[9-11].  The  d^ree  to  which  the  levels  of  these  other 
complexes  alter  the  physical  distribution  and  activity  of 
Myc  at  E-box  versus  non-E-box  promoters  remains  to  be 
rigorously  evaluated.  The  si)ecialized  r^ulation  of  these 


other  E-box-binding  proteins  (such  as  Mondo)  have  the 
capacity  to  tune  the  Myc  response  toward  or  away  fi'om 
cellular  processes,  such  as  metabolism  [59]  or  circadian 
rhythms  via  coupling  with  proteins  such  as  CLOCK  [60]. 

Almost  certainly,  steady-state  patterns  of  gene  expres¬ 
sion  observed  after  Myc  expression  do  not  only  reflect  the 
direct  action  of  Myc,  but  are  also  due  to  secondary  effects  of 
proteins  or  RNAs  encoded  by  Myc  target  genes  that  feed 
back  on  gene  expression  [38] .  In  a  rapidly  growing  cell,  such 
as  a  B  cell  becoming  activated  fi:um  the  naive  state,  increas¬ 
ing  s5mthesis  of  macromolecules  is  required  to  satisfy  a 
growing  demand  for  the  transcription  machinery  and  other 
components  required  to  activate  genes.  Examples  for  such 
indirect  effects  are  the  induction  of  the  (jcn5  histone  acet- 
ylase  that  mediates  a  global  opening  of  chromatin  upon  Myc 
expression  and  enhanced  expression  of  phosphoribosyl-py- 
rophosphate  synthetase  2  (PRPS2),  which  promotes  the 
enhanced  nucleotide  biosynthesis  of  Myc-driven  tiunors 
[18,26,61].  At  the  same  time,  rapid  growth  and  enhanced 
global  transcription  may  induce  the  redeployment  of  limit¬ 
ing  components  of  the  gene  expression  machinery  to  more  or 
less  active  transcription  start  sites.  Under  conditions  of 
limiting  gene  exjjression  components,  increasing  Myc  at 
strongly  expressed  promoters  might  recruit  transcripjtion 


Figure  3.  Schematic  showing  the  dynamic  occupancy  of  E-box  sites  by  Myc  and  other  transcriptional  regulators.  Multiple  basic  helix-loop-helix  (bHLH)  transcription  factor 
pairings  compete  for  binding  at  E-box  sites.  The  net  positive  or  negative  transcriptional  consequence  of  regulator  binding  at  E-boxes  is  depicted.  For  definitions  of 
abbreviations,  please  see  the  main  text. 
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Rgura  4.  Schematic  of  Myc  transfer  functions.  Transfer  functions  are  shown  providing  the  relation  between  expression  changes  at  genes  as  a  function  of  various 
parameters  as  measured  either  by  absolute  or  relative  change  in  expression.  (A|  A  linear  amplifier  transfer  function  where  the  change  in  expression  in  high  versus  low  Myc 
conditions  is  constant  across  all  transcriptionally  active  genes.  IB|  A  nonlinear  amplifier  transfer  function  where  the  change  in  expression  in  high  versus  low  Myc  conditions 
is  pnoport'nnal  to  the  exponent  of  the  expression  in  low  Myc  conditions.  IC|  A  Myc/Mizl  repression  transfer  function  where  the  change  in  expression  is  proportional  to  the 
ratio  of  AMyc  binding  divkied  by  Mizl  binding.  |D|  A  combined  Mizi-shaped  nonlinear  amplifier  model  that  combines  the  transfer  functions  in  (B)and  (C).  A  2D  heatmap  is 
provided  showing  the  relative  change  in  expression  as  a  function  of  both  the  ratio  of  AMyc  binding  divided  by  Mizl  binding  (y-axis)  and  the  initial  expression  in  low  Myc 
conditions  (x-axis). 
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machinery  otherwise  assigned  to  weakly  transcribed  genes. 
Such  diversion  would  yield  the  upregulation  of  some  genes 
with  the  downregulation  of  others.  This  repression  by  redis¬ 
tribution  of  the  transcriptional  machinery  has  been  recently 
observed  during  nuclear  factor  (NF)-kB  activation  in  endo¬ 
thelial  cells  [62]. 

In  principle,  other  mechanisms  may  act  to  counterbal¬ 
ance  transcriptome  amplification.  For  example,  in  yeast, 
rates  of  RNA  synthesis  and  RNA  degradation  are  coupled 
[63] .  Mutations  and  strategies  that  alter  the  levels  of  one 
elicit  compensatory  changes  in  the  other.  Additional  am¬ 
plification  counteracting  pressures  may  arise  from  rate 
limiting  levels  of  translation  initiation  in  tumors.  Notably, 
disruption  of  the  translation  initiation  factor  eukaryotic 
translation  initiation  factor  4  (eIF4F)  is  S5mthetic  lethal  in 
tumors  with  deregulated  Myc,  suggesting  that  the  high 
levels  of  mRNA  output  need  to  be  supported  by  high  levels 
of  protein  expression,  perhaps  exposing  an  orthogonal 
vulnerability  arising  due  to  amplification  [64].  For  RNA 
levels  to  rise  globally  in  response  to  MYC  deregulation, 
these  compensatory  countermeasures  may  need  to  be 
defeated  or  overwhelmed. 

Concluding  remarks 

Deregulated  in  tumors,  Myc  acts  as  a  universal  regulator  of 
the  pre-existing  gene  expression  program  of  the  cell.  Wheth¬ 
er  amplification  alone  provides  a  sufficient  oncogenic  im¬ 
pulse  or  whether  Myc-dependent  repressive  mechanisms 
are  required  for  Myc-driven  tumorigenesis  warrants  further 
investigation.  Simply,  as  a  universal  regulator,  mutations  in 
other  oncogenes  and  tumor  suppressors  must  first  specify 
particular  genes  to  be  turned  on  or  off,  and  Myc  then 
develops  the  new  gene  expression  program.  This  model  is 
also  consistent  with  data  in  iPS  reprogramming,  where  Myc 
overexpression  increases  the  efficiency  of  an  embryonic  stem 
cell  specif3dng  Octamer  4/SRY  (sex  determining  region  Y)- 
box  2/Kruppel-like  factor  4  (Oct4/Sox2/Klf4)  transcription 
factors  [65] .  Alternatively,  the  pathways  for  growth,  prolif¬ 
eration,  apoptosis,  and  neoplastic  metabolism  do  not  scale 
linearly  (i.e.,  increasing  the  components  of  these  pathways 
drives  them  across  thresholds  that  gate  their  activity  as  a 
step  function).  In  this  situation,  malignant  behavior  would 
manifest  as  an  emergent  property  explainable  by  the  control 
theory  (Box  1)  and  architecture  of  subsystems  and  networks 
whose  components  are  upregulated  by  Myc. 

Although  predominantly  an  amplifier  of  transcription, 
context-specific  conditions  and  constraints  at  individual 
genes  modulate  Myc  transcriptional  output  and  response. 
To  elucidate  the  molecular  mechanism  of  Myc-driven  tran¬ 
scriptome  amplification,  a  better  description  of  its  transfer 
function  (the  function  describing  how  the  Myc  transcrip¬ 
tional  effect  depends  on  initial  parameters  at  any  gene;  i.e., 
the  input-output  curve)  will  be  required  to  account  for  its 
physiological  and  pathological  activities.  If  proportional 
and  linear,  for  example,  all  Myc  targets  will  be  upregulated 
equivalently,  and  any  changes  in  the  relative  expression  of 
particular  genes  would  demand  parallel  changes  in  the 
levels  of  the  gene-specific  factors  that  differentially  up-  or 
downregulate  these  targets  (Figure  4A).  However,  if  the 
transfer  function  is  nonlinear,  as  proposed  in  [66],  then 
different  zones  across  the  expression  spectrum  will  be 


Box  1.  Control  theory 

Control  theory  is  a  dynamic  analysis  of  network  behavior  as  modeled 
by  the  effect  of  inputs  and  their  consequent  feedback  on  system 
output,  in  engineering,  this  can  describe  how  different  components 
affecttheoutput  of  a  system.  For  instance,  a  compressor  either  boosts 
or  dampens  signal  depending  on  initial  input  signal.  A  reverb  module 
would  apply  a  delay  to  an  input  signal  before  entering  it  back  into  the 
system,  creating  a  recursive  feedback  loop.  An  amplifier  increases 
signal  by  an  amount  relative  to  the  initial  input  signal.  In  each  of  these 
examples,  the  transformations  are  surmised  through  the  transfer 
function:  a  mathematical  representation  of  the  relation  between  input 
and  output.  Thus,  a  control  theory  analysis  of  Myc  examines  the 
consequence  of  adding  Myc  to  the  transcriptional  system.  The  Myc 
transfer  function  describes  the  transcriptional  output  of  Myc  activity 
as  a  function  of  the  initial  transcription  parameters. 

disproportionately  regulated.  If  highly  expressed  genes 
are  preferentially  amplified,  then  a  Myc  signature  may 
emerge  comprising  genes  that  are  highly  expressed  across 
many  tissues  (Figure  4B).  Of  course,  such  amplification 
cannot  continue  indefinitely  and  unopposed;  supraphysio- 
logical  levels  of  Myc  are  expected  to  elicit  either  compen¬ 
satory  measures  or  cell  death  (Figure  4C).  Increasing 
amplification  of  all  genes  demands  that  the  cell  can  mar¬ 
shal  sufficient  resources  to  augment  net  transcription.  In 
the  presence  of  limiting  materials,  the  competition 
between  promoters  may  modify  the  transfer  function  from 
one  that  supports  only  universal  amplification  to  one  that 
allows  for  repression  at  those  promoters  deprived  of 
limiting  factors  that  have  been  deployed  to  support  the 
expression  of  other  targets.  In  this  scheme,  repression  (e.g., 
Miz  1-dependent  repression)  reflects  feedback  control  that 
responds  to  the  unbalanced  and  excessive  Myc-enforced 
expression  of  other  targets  (Figure  4D). 

Universal  activation  by  Myc  does  not  preclude  the  func¬ 
tional  emergence  of  preferred  and,  hence,  specific  targets. 
Moreover,  the  preferential  up-  or  downregulation  of  specif¬ 
ic  genes  that  occurs  as  tumors  achieve  homeostasis  in  high 
Myc  conditions  may  underlie  the  ability  of  Myc  to  act  as  one 
of  the  rare  universal  oncogenes  capable  of  promoting 
tumorigenesis  in  a  wide  variety  of  tissues.  Only  when 
experiments  challenge  models  that  make  quantitative, 
precise,  and  testable  predictions  will  a  durable  consensus 
on  the  role  of  Myc  emerge.  Beyond  having  heuristic  value 
as  an  organizing  framework  for  diverse  observations,  the 
universal  nonlinear  amplifier  model  may  yet  account  for 
quantitative  differences  in  gene  expression  between  cell 
states  in  health  and  disease. 
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